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CN I ABSTRACT 

i— I ■ I calculate the physical properties of 32 transiting extrasolar planet and brown-dwarf systems 

from existing photometric observations and measured spectroscopic parameters. The systems 

studied include fifteen observed by the CoRoT satellite, ten by Kepler and five by the Deep 

^vO • Impact spacecraft. Inclusion of the objects studied in previous papers leads to a sample of 58 

transiting systems with homogeneously measured properties. The Kepler data include obser- 
^ , vations from Quarter 2, and my analyses of several of the systems are the first to be based on 

short-cadence data from this satellite. 

The light curves are modelled using the JKTEBOP code, with attention paid to the treat- 
SZ . ment of limb darkening, contaminating light, orbital eccentricity, correlated noise, and nu- 

Q-T merical integration over long exposure times. The physical properties are derived from the 

q \ light curve parameters, spectroscopic characteristics of the host star, and constraints from five 

Vh ■ sets of theoretical stellar model predictions. An alternative approach using a calibration from 

r/} ' eclipsing binary star systems is explored and found to give comparable results whilst imposing 

d . a much smaller computational burden. 

My results are in good agreement with published properties for most of the transiting 

, \ systems, but discrepancies are identified for CoRoT-5, CoRoT-8, CoRoT-13, Kepler-5 and 

^. , Kepler-7. Many of the errorbars quoted in the literature are underestimated. Refined orbital 

iy~j ' ephemerides are given for CoRoT-8 and for the Kepler planets. Asteroseismic constraints 

(T) \ on the density of the host stars are in good agreement with the photometric equivalents for 

(N| ■ HD 17156 and TrES-2, but not for HAT-P-7 and HAT-P-11. 

t-H ' Complete error budgets are generated for each transiting system, allowing identification 

r — | , of the observations best-suited to improve measurements of their physical properties. Whilst 

f*^ ■ most systems would benefit from further photometry and spectroscopy, HD 17156, HD 80606, 

j— | ' HAT-P-7 and TrES-2 are now extremely well characterised. HAT-P- 1 1 is an exceptional candi- 

t— I . date for studying starspots. The orbital ephemerides of some transiting systems are becoming 

LJ uncertain and they should be re-observed in the near future. 

. £h ! The primary results from the current work and from previous papers in the se- 

S^ . ries have been placed in an online catalogue, from where they can be obtained 

Vh ' in a range of formats for reference and further study. TEPCat is available at 

|http : //www. astro. keele.ac .uk/~ jkt/tepcat/| 

Key words: stars: planetary systems — stars: binaries: eclipsing — stars: binaries: spectro- 
scopic — stars: fundamental parameters 



1 INTRODUCTION selection effects mean that most known TEPs orbit very close to 

m their parent star and thus have highly irradiated atmospheres whose 

Of the 550 planets known to orbit stars other than our Surrl the , . , . , ..,.,., .. 

physical properties can be scrutinised using high-quality astronom- 
transiting systems are the most interesting ones. Transiting extra- , , 

, ™™, . , , . , „ ? ical observations, 

solar planets (TEPs) are the only planets outside our Solar sys- 
tem whose masses and radii, and thus surface gravities and mean One hindrance to the study of TEPs is that measurement of 
densities, can be measured to reasonable precision. Observational their physical properties requires not only transit light curves and 

radial velocity measurements of the parent stars, but also some sort 

of additional constraint. This is ordinarily obtained by forcing the 

* E-mail: jkt@astro.keele.ac.uk properties of the star to match predictions from theoretical stellar 

1 www . exoplanet . eu models, guided by an effective temperature (T c ft) and metal abun- 
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dance ([§]) obtained from spectral analysis. The dependence on 
stellar theory leads to systematic errors which can be sizeable for 
some of the measured quantities, and also allows inhomogeneities 
to occur between studies which use different theoretical predictions 
or apply the constraint in a different way. This in turn compromises 
statistical studies of transiting planets. 

For these reasons I am measuring the properties of the known 
TEPs using strictly homogeneous methods. Paper I ( iSouthworthl 
2008) discussed the methodology used to model the transit light 
curves, paying particular attention to error analysis and the treat- 
ment of limb darkening, and applied them t o the 14 systems with 
good observational data at the time. Paper II (Southworth 2009) ex- 
plored the application of constraints from seven different sets of 
theoretical model predictions and alternatively an empirical mass- 
radius relation obtained from eclipsing binaries. Three of the sets 
of theoretical predictions were selected to aid in the determination 
of the physical properties of the same 14 systems, with detailed er- 
ror bud gets including ra ndom and systematic contributions. In Pa- 
per III dSouthworthluOlOh I extended the number of systems to 30 
and the number of theoretical predictions used to five sets, improv- 
ing both the statistical weight of the ensemble and the precision of 
the systematic errorbars. 

In the current work I enlarge the number of TEPs with ho- 
mogeneous properties to 58, concentrating on those which hav e 
been observed by the space missions Kepler (Borucki et al. 2010b), 
CoRoT JBaglin et al.l2006t) and EPOCH dChristiansen et al.l2 009). 
In Paper II I outlined the concept of applying an observational 
mass-radius relation instead of using constraints from stellar the- 
ory, resulting in totally empirical measurements of the properties of 
TEP systems. This approach was not very successful (see Paper III) 
because the mass-radius relation did not allow for stellar evolu- 
tion and had to be calibrated on low-mass eclipsing binary systems 
whose properties (primarily radius) are affected by magnetic activ- 
ity arising from comparatively f ast rotation. In the c urrent work I 
follow the alternative approach of lEnoch et alj d2010f) which relates 
the radius of the star to its density, T e ff and [§] • This technique 
is not purely empirical, as it incorporates parameters derived using 
spectral synthesis techniques for both the TEP hosts and the cali- 
brating sample, but returns results in much better agreement with 
the default method using theoretical predictions. Finally, I explore 
the opportunities for further study of the systems studied in this 
work. 

For a small number of TEPs it is possible to either partially 
or totally avoid systematic errors: the study of transit timing vari- 
ations (T TVs) allow the masses o f some TEPs to be co nstrained 
directly JHolman & Murray! 120051 : iLissauer et ail 1201 lh : and ra- 
dial velocitie s of the planet HD 209458 b have been measured by 
ISnellen et alj ( 1201 Oaf) from infrared absorption lines, allowing the 
system properties to be calculated in an identical way to double- 
lined eclipsing binary star systems. These methods are, however, 
not applicable to the predominant population of TEPs which show 
no detectable TTVs and are not amenable to direct velocity mea- 
surements. 



2 ANALYSIS OF THE LIGHT CURVES 

I have modelled the light curves of each TEP using the methods 
espoused in Paper I. In short, the JKTEB OBJ code {Southworth et alj 



2 JKTEBOP is written in FORTRAN77 and the source code is available at 

|http : //www . a stro . kee le .ac .uk/^j k t /codes / jktebop .html 



l2004aU bl) is used to model the available transit light curves. The 
components are approximated as biaxial spheroids whose shapes 
are governed by the mass ratio, q. The results in this work are all 
extremely insensitive to the values adopted for q. 

The main parameters of a JKTEBOP fit are the orbital inclina- 
tion, i, and the fractional radii of the star and planet, ta and r^. 
These are denned as 



r A = 



Ra 



i-h 



Rb 

a 



(1) 



where Ra and Rh are the volume-equivalent stellar and planetary 
radii and a is the orbital semimajor axis. In JKTEBOP the fractional 
radii are re-parameterised as their sum and ratio: 



rA + r h 
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as these are less strongly correlated. In general the orbital period, 
Porb, is taken from the literature and the time of transit midpoint, 
To, is included as a fitted parameter. 

Each light curve is fitted with a number of different ap- 
proaches to limb darkening (LP ), la errorbars are obtained usin g 
1000 Monte Carlo simulations dSoufhworth et al.ll2004cl , l2005bl) . 
Errorbars are also calcul ated using a resid ual permutation (or 
"prayer bead") algorithm ( Jenkins et alj|2002i) which accounts for 
correlated observational noise, and the largest of the two alterna- 
tives is adopted for each parameter. 

The LD of the star has an important influence on transit light 
curves. For each light curve, solutions are obtained using five differ- 
ent LD laws, each with three different approaches to the limb dark- 
ening coefficients (LDCs): (1) both fixed (hereafter 'LD-fixed'); (2) 
the linear one (ma) fitted and the nonlinear one (va) fixed but per- 
turbed by ±0.05 in the error analysis simulations ('LD-fit/fix'); and 
(3) both coefficients fitted ('LD-fitted'). Initial or fixed values for 
the LDCs are bilinearly interpolated in T a g and logo within the 
tabulated theoretical predictions included in the JKTLlJj code. 

Once the best of the three alternatives (LD-fixed, LD-fit/fix, 
LD-fitted) is identified, the combined solution for that option is 
calculated by taking the mean of the solutions for the four two- 
coefficient LD laws and by taking the largest errorbar from these 
solutions plus a contribution to account for the scatter in the pa- 
rameter values. In most cases the LD-fit/fix solutions turn out to 
be the best compromise between severing the dependence on the- 
oretical calculations and trying to fit too many parameters to the 
data. 

Some TEPs have a non-circular orbit which changes the dura- 
tion of the trans it but has a negligible effect on the light curve shape 
dKippindl2008r) . I account for this by adding published constraints 
on orbital eccentricity (e) and periastron longitude (lu), with the pa- 
rameter combinations e cos a; and e sin ui when possible, us ing the 
approach outlined in Paper III and lSouthworth et alj (12009c). 

An additional annoyance for some systems is light from a 
nearb y star contaminating the photometry (e.g. iDaemgen et alj 
2009). This 'third light' makes the transit shallower but cannot be 
fitted for in the light curve due to strong degeneracy with ta, r^ and 
i (see Paper III). When third light is known to exist it i s acco unted 
for using the method put forward by S outhworth et alj (2010). 
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Figure 1. Parameters of light curve solutions of an undersampled light 
curve, closely resembling that of Kepler-6, with different numbers of nu- 
merical integration points used in the solution. Unbroken lines show the 
correct parameter values and dotted lines show the ± 1 % intervals for these 
parameters. 
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Figure 2. Plot of a model light curve closely resembling that of Kepler-6, 
but covering 32 instead of 12 orbital periods (blue points). This is compared 
to the best fit obtained without using numerical integration (red line). The 
residuals of the fit are plotted at the base of the figure offset from zero. 



2.1 Numerical integration of light curves 

In some light curves the sampling rate is a significant fraction of 
the transit duration, leading to a smearing of the transit shape. If 
left uncorrected this could cause errors of up to 30% (worst-case 



3 JKTLD is written in FORTRAN77 and is the source code is available at 

|http: //www .astro .keele .ac . uk/~ jkt / codes/ jktld.html 



scenario) in the physical properties of TEPs. The prime candidates 
for this problem are the Kepler planets, whose long-cadence photo- 
metric points consist of summations of 270 consecutive datapoints 
leading to an overall s ampling rate of one datum every 29.4244 min 
jjenkins et alj|2010bh . Some of the CoRoT satellite data are also 
affected, as the standard cadence for this instrument is 512 s. Most 
CoRoT TEPs also have short-cadence sampling with a rate of 32 s; 
these data do not suffer from temporal undersampling. 

For the current work I have modified JKTEBOP to optionally 
perform numerical integration to cope with the Kepler and CoRoT 
data. The approach is to calculate a given number of model dat- 
apoints (iVi n t) evenly spread over a given time interval, and sum 
them to create an integrated datapoint which can be directly com- 
pared with observations. The next question is: how fine a time sam- 
pling is necessary? I created a dataset very similar to that of the 
TEP Kepler-6 by generating a model light curve (r\ + T~b — 0.16, 
k — 0.10, i = 86.5°, quadratic LD), extending it over 12 orbital 
cycles, and summing it into 29.4 min bins. This was then fitted with 
/Vint varied from 1 (i.e. no numerical integration) to 15 (equivalent 
to a 1.96 min sampling rate). The resulting values of va + ft>, k 
and i are plotted in Fig.[T]and show that ta + Tb is more affected 
than k and i. Using 7v"i nt = 10 means we incur an error of only 
0.1% in r\ + r\,', using iVint = 5 would engender a 1% error. 
TEPs with shorter orbital periods (i.e. quicker transits) or higher i 
(sharper partial phases) will be more strongly affected. 

Fig-El shows the Kepler-6-like model light curve extended to 
cover 32 orbital cycles and summed into 29.4 min bins. This has 
been fitted by JKTEBOP but without performing numerical integra- 
tion, in order to demonstrate the effect of neglecting the undersam- 
pling. The fitted model is unable to correctly reproduce the syn- 
thetic data during the partial phases of the transit, the discrepancy 
being worst at the first and last contact points where the light curve 
derivative is of greatest magnitude. This suggests that it would be 
possible to fit for the amount of numerical integration needed for a 
high-quality dataset, although it is very unlikely that such an option 
will ever be useful. The amount of numerical integration needed 
will also be quite sensitive to LD, particularly around the limb of 
the star. 



3 CALCULATION OF PHYSICAL PROPERTIES 

3.1 Via stellar models 

Analysis of a transit light curve gives the quantities P or b> To, ta, rb 
and i. Radial velocity (RV) measurements of the parent star yield its 
orbital velocity amplitude, Ka- Measuring the physical properties 
of the system requires an additional constraint, which is normally 
taken from theoretical stellar evolutionary models. The observed 
stellar effective temperature, T e e and metal abundance, [W] , are 
useful to guide this process, which is discussed in detail in Paper II 
and Paper III. 

I use the velocity amplitude of the planet (Kb) to control the 
solution process. A starting value is guessed, and is combined with 
the measured P or b, ?"a, J"b, i and Ka to obtain the physical proper- 
ties of the system using standard formulae (e.g. Hilditch 2001). The 
resulting calculated stellar mass (Ma) and measured [^1 are used 
to obtain the expected stellar radius (Ra) and T a g by interpolation 
in a set of tabulated theoretical predictions. Kb is then iteratively 
refined to obtain the best agreement with the calculated Ra and 
observed T c ff by minimising the figure of merit: 
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This process is performed for a range of ages from the zero-age to 
the terminal-age main sequence (curtailed at a maximum of 20 Gyr) 
to find the overall best fit. The code which performs this step (JK- 
TABSDIM) has been profiled in order to improve its speed, which 
has allowed the step size in age to be decreased from 0. 1 Gyr to 
0.01 Gyr for the current work. 

The uncertainties on the input para meters to JKTABSDIM ar e 
propagated using a perturbation analysis dSouthworth et al.l2 005a). 
resulting in a complete error budget for every output parameter. 
This allows identification of which type of observations would be 
best to improve our understanding of each TEP (see Sect.[9). 

Apart from the random errors, which are calculated using the 
propagation analysis, systematic errors arise from the use of theo- 
retical stellar models. These can be estimated by running solutions 
with a range of different model predictions. As in Paper III I use five 
set s of model predictions: Claret (Claret 2004, 2005, 2006, 2007), 
F 2 dDemarque et alj|2004h . Teramo jPietrinferni et alj[2 004). VRSS 
(VandenBerg et al .1 120061) and DSEP JDotter et al.ll2008nT 7KTABS- 
DIM is run with each of these five sets, and the final results are 
taken to be the unweighted mean of the individual values for each 
output quantity. The statistical error is taken to be the largest of 
the individual uncertainties from the perturbation analysis, and the 
systematic error to be the standard deviation of the values from 
each of the model sets. The final results therefore rest evenly on the 
predictions of all five model sets. Since Paper III I have obtained 
additional tabulations for the Claret models, for fractional metal 
abundances of Z = 0.005 and Z = 0.015, to allow for [§] val- 
ues down to —0.60. I have also identified and fixed a mistake in 
my reformatting of the VRSS models which caused the wrong Ma 
values to be used for a small number of tabulations. 

The final results from five runs of JKTABSDIM with different 
theoretical model sets is the following parameters with statistical 
and systematic errorbars: the mass, radius, surface gravity and den- 
sity of the star (Ma, Ra, log qa, Pa) and of the planet (Mb, Rb, 
<7b, Pb)- In additional to this I calculate a surrogate for the equilib- 
rium temperature for the planet: 

and also its lSafronovl d 19721) number: 

I 1 = I — I I ■ '' 
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Three quantities are indepe ndent 
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dSouthworth eTail l2007bh , p A dSeager & Mallen-Omelasl 120031) 
and T e ' q (Paper III). 

3.2 Via eclipsing binary relations 

In Paper II I found a way to bypass the use of stellar models 
entirely, by defining an empirical mass-radius relation based on 
well-studied and unevolved detached eclipsing binary star systems 
(dEBs). The chief advantages of this approach were tractability and 
the avoidence of a dependence on stellar theory. The primary disad- 
vantage was that the results were much inferior to those calculated 
using stellar models. This arose because the known well-studied 
dEBs tend to h ave substantially larger radii than predicted by the 
models (see e.g. Lopez-Morales 2007; Ribas et al. 2008), so a more 
massive star was needed to reproduce the density obtained from 
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Figure 3. Plot of density versus T c g for TEP host stars compared to well- 
studied dEBs. Filled (red) circles represent data from Paper III and the cur- 
rent work and open circles show data from the literature. Plusses and crosses 
(blue) denote the primary and secondary components of dEBs, respectively. 



the light curve solutions. Stellar evolution was also not allowed for, 
so the systems with more evolved stars were found to be rather 
more massive than they actually are. These problems make a sim- 
ple mass-radius relation untenable. 

An alternative approa ch would be to use the empirical rela- 
tions for dEBs defined by iTorres et alj d2010f) . which give log M 
and logR in terms of logT of f, logg and \W\. These account 
for evolution and also for metal abundance, and have a modest 
scatter o f a = 0.02 7 in lo gM and a = 0.014 in logi?. As 
noted by lEnoch et al] d20100 . a better approach for TEPs would 
be to replace logy with logp as the former quantity is rather 
tricky to derive from spectra of solar-like stars, whereas the lat- 
ter quantity is almost d irectly measura ble from t ransit light curves 
dSeager & Mallen-Ornelas] |2003|) . lEnoch et al] d201Cj|) found that 
this modified approach yields lower scatters of a = 0.023 in log M 
mid a = 0.009 in log R. 

There are se veral possible criticisms of the implementation 
of this method by lEnoch et al] ([2010). Firstly, they included only 
those dEBs with metallici ty measurements ( 19 out of the 95 sys- 
tems in the compilation bv lTorres et alj|2010l) so could suffer from 
small-number statistics. Also, many of the TEP host stars are of 
mass <1.3 Mq where metal abundance has a much smaller effect 
on the stellar propertieQ Secondly, they included component stars 
of dEBs with masses up to 15 Mq , which are of uncertain value for 
studying the currently know TEP hosts. Finally, a plot of density 
versus T e s does not allow one to conclude that the dEBs are good 
calibrators of TEP host stars: the former might be systematically 
less dense than the latter (see Fig.[5J. This leads to concerns that 
the relations have to b e extrapolated to TE P hosts rather than inter- 
polated. However, the lEnoch et a l. (2010) approach yields results 
which are quick to calculate and are in good agreement with stellar 
theory, so is worthy of further investigation. 

I have therefore obtained my own calibration of log 10 M and 



4 The value of 1.3 Mq below which metal abundance is comparatively 
unimportant was obtained by plotting VRSS zero-age main sequence 
isochrones for Z values ranging from 0.01 to 0.05 
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Table 1. Coefficients of the equations for log 10 M and log 10 R derived using eclipsing binary star systems. Extra significant figures are provided to guard 
against round-off errors. 



Calibration parameter 


log 10 M 


logic M 


logio R 


log 10 R 


Valid mass range ( Mg) 


0.2 to 27.0 


0.2 to 3.0 


0.2 to 27.0 


0.2 to 3.0 


Cl 


0.01092 ±0.00176 


0.01384 ±0.00263 


0.003759 ± 0.000933 


0.007055 ± 0.000930 


CI 


1.0826 ±0.0178 


1.0569 ±0.0268 


0.36325 ±0.00711 


0.37796 ±0.00729 


C3 


0.4028 ± 0.0259 


0.360 ±0.107 


0.1317 ±0.0102 




C4 


-0.15139 ±0.00536 


-0.16236 ±0.00936 


-0.38296 ±0.00161 


-0.37880 ±0.00175 


C5 


-0.0008 ±0.00159 


-0.01405 ±0.00589 






c% 


0.1803 ±0.0105 


0.1755 ±0.0102 


0.06020 ±0.00417 


0.06004 ± 0.00422 


Scatter (dex) 


0.0286 


0.0268 


0.00953 


0.00907 



log 10 R in terms of log 10 T e g , log 10 p and y^\ ■ The calibration 
sample is taken from DEBCajjand includes all stars but onqjwith 
masses up to 3Mq. Objects without a metallicity measurement 
were assigned [^1 = ± 100. The final sample contains 90 dEBs 
(180 stars), and benefits in partic ular from the six new low-mass 
dEBs studied bv lKraus et aUd201lh . 
The adopted equations are: 

-.2 



Table 2. Corrected planetary densities of the TEPs studied in Paper III. 



log 10 M = ci + c 2 log 10 X + c 3 (log 10 X) 

+C 4 lOg 10 p + C 5 (log 10 p) 2 + C 6 



logi R = ci+c 2 log 10 X + c 3 (log 10 X) 

+C 4 logi p + C 5 (log 10 pf + C 6 



(6) 



(7) 



where X = log 10 (T c ff /T D s,rs>), Toff,© = 5781 K, and mass, radius 
and density are given in solar units. The coefficients of the fit, a, 
were found using a downhill simplex algorithm and are given in Ta- 
bleQ] The equation adopted for the current paper is for log 10 R with 
a range of validity of 0.2-3.0 M©. The scatte r around the calibra - 
tion is only 0.009 dex, in good agreement with lEnoch et al.l (2010). 
The coefficients are much more precise because of the larger cal- 
ibration sample, but are not directly comparable because of a dif- 
ferent choice of normalisation parameter for T a g . The calibration 
has been implemented into the JKTABSDIM code, and its scatter 
is propagated through into the final results using the perturbation 
method. 



3.3 Constants and units 

In previous papers in this serie^j the densities of TEPs were given 
relative to that of Jupiter. The density of Jupiter was calculated 
incorrectly, using the equatorial radius rather than the volume- 
equivalent radius. The former is larger by 2.26% due to the oblate- 
ness of the planet, leading to a scale in which the density of the 
planet Jupiter was counterintuitively 1.0694 pj up . Starting with the 
present work I account for this effect, leading to planetary densities 
which are lower by 6.94%. No other quantities are affected by this 
oversight. In Table|2]l give the corrected densities of those planets 
which were studied in previous papers (see Paper III). It is likely 

5 Catalogue of well-studied detached eclipsing binary star systems: 
|http : //www .astro .keele .ac.uk/^jkt/debcat/| 



System 



Pb ( PJup) 



GJ436 

HAT-P-1 

HAT-P-2 

HD 149026 

HD 189733 

HD 209458 

OGLE-TR-10 

OGLE-TR-56 

OGLE-TR-111 

OGLE-TR-113 

OGLE-TR-132 

OGLE-TR-182 

OGLE-TR-21 1 

OGLE-TR-L9 

TrES-1 

TrES-2 

TrES-3 

TrES-4 

WASP-1 

WASP-2 

WASP-3 

WASP-4 

WASP-5 

WASP- 10 

WASP- 18 

XO-1 

XO-2 

XO-3 

XO-4 

XO-5 



141+018 
-0.02 
n 971 +0.025 
u - zll -0.002 

4 8+1-4 
4 '°-0.0 

-, c-7+0.67 +0.01 
1 -°' -0.53 -0.01 

n 7nfi+ - 062 + 0008 

U.IUO _ 062 -0.008 

254+ ' 004 
u.zoi_ 002 

n 19 c;+0.036 
u - lzd -0 .000 

u -' u _0.32 

0-85+S:oi 
u -°°-o.oo 

311 +0.104 
u - oll -0.002 

n q/io +0-109 +0.001 

u.o<iO _ 124 -0.001 

96+ ' 35 
u - ao -0.01 
5Sfi+0052 
u.ojo_ ()09 

U.t)84_ 0006 

n SO/1 +0-047 +0.006 
u.out_ 053 _ .004 
n 1QS+0.037 

u.id«_ 0001 

n OAF. +0-055 +0.001 
U.Atv _ .033 -0.001 

70+ ' 13 
u -' u -o.oi 

63+ 010 
u - DC> -o.oo 

n zifiq+o- 014 +0.005 

u.ioo _ .017 -0.006 

QQO + 0.13 

u - 3 °-o.oi 

o ^q+0.36 
z -^°-0.05 

6 21+ ' 84 

D - Zi -0.05 

4Q2+0-059 

"•*"-() .005 

KQQ +0.111 +0.012 

U.JOi _o.o67 -0.010 

5 69+ ' 63 

°-" 3 -0.06 
n fifi +0.12 +0.01 
u - uu -0.37 -0.00 
79+ 017 

u -' 3 -o.oo 



ASAS J0528+0338 was dropped from the sample as it is a pre-main- 

sequence system. 

7 Paper II and Paper III, and also lSouthworth et all b009jblcll2oTol|201lh . 



that this error has occurred several times before in the literature, 
and in these cases pb should be divided by 1.0694. In the Appendix 
I make no attempt to perform this correction to published values as 
it is not clear which of them are affected. 

This seems an appropriate point to specify the values of all 
constants and adopted quantities used in the present series of pa- 
pers. The measured physical properties of the star are given in solar 
units for mass, radius and mean density, and in c.g.s. for log g. For 
the planet, mass, radius and density are given in Jupiter units, and 
surface gravity in ms~ 2 . Table[3]gives all constants used, plus the 
ensuing values of some quantities of interest. Since this paper was 
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Table 3. Physical constants and adopted quantities used in the current work. 
The lower part of the Table gives the value of several quantities resulting 
from the adopted physical constants. 

References: (1) US National Institute of Standards and Techno logy (2006 
constants); (2) IERS Conve ntions 2010 IPetit &Luzumll20irj|); (31 2010 
Astronomical Alma nac; (41 iBrown & Chri stensen-Dalsgaard (1998); (5) 
iBahcall et al] {\99JY . (6) JPL DE405 Ephemerides; (7) NASA NSSDC 
Jupiter fact sheet. 



Parameter 



Value 



Unit 



Ref. 



Stefan-Boltzmann constant 


5.67040 X 10" 8 


Wm~ 2 K" 4 


1 


G 


6.67428 x 10~ n 


m 3 kg -1 s _1 


1 


GM e 


1.32712421 x 10 20 


m 3 s -2 


2 


AU 


1.49597871 x 10 11 


m 


3 


Rq 


6.95508 x 10 9 


m 


4 


Lq 


3.844 x 10 26 


W 


5 


i?j up (equatorial) 


7.1492 x 10 8 


m 


3 


M /M Jup 


1.0473486 x 10 3 




6 


Volume of Jupiter 


1.43128 x 10 24 


m 3 


7 


M Q 


1.98842 x 10 30 


kg 




logfl© 


4.43831 


(c.g.s.) 




P0 


1410.95 


kgm -3 




Mj U p 


1.89852 x 10 27 


kg 




Sjup (equatorial) 


24.7916 


ms -2 




Pjup 


1326.45 


kgm -3 





submitted, lHarmanec & Prsal u201 lh have proposed a standardisa- 
tion of the physical constants used in astronomy. They propose that 
a nominal solar mass, radius and luminosity are defined as exact 
quantities for use by all researchers. The values they propose are in 
full agreement with those used in the current work. 



4 LIGHT CURVE MORPHOLOGY 

The literature contains a proliferation of terms for the different 
types of eclipses in TEP systems. The correct terminology has 
been establish ed for many years for eclipsing binary systems (e.g. 
iHilditchl 1200 lh - and a TEP is simply a special case of an eclips- 
ing binary. A 'transit' is when a smaller object passes in front of 
a larger object, for example a planet in front of a star. An 'oc- 
cultation' is when a smaller object passes behind a larger object. 
Transits and occultations are two types of eclipses, and the third 
type is a 'partial eclipse' . Fig.|4]is a schematic representation of the 
situation. 

The definition of 'primary eclipse' is when the hotter object 
is behind the cooler object, and 'secondary eclipse' denotes to the 
reverse situation. This ordinarily means that the primary eclipse is 
deeper than the secondary, although exceptions to the rule are pos- 
sible in the case of eccentric orbits (where different surface areas 
are covered at the two types of eclipse). For almost all TEPq pri- 
mary eclipses are transits and secondary eclipses are occultations. 

Eclipsing systems can have only one transit per orbijj, so ref- 
erences to a 'primary transit' or 'anti-transit' are misleading. The 
phrase 'secondary transit' is incorrect and should not be used. 




Secondary eclipse 

V 



Occultotion 



0.0 



Transit 
0.2 



0.4 0.6 

Orbital phase 



0.8 



Figure 4. Representation of a transit light curve with the main features in- 
dicated and labelled. 



5 DATA ACQUISITION 

The C0R0T data used here are the N2 public version obtained from 
the public archiv£3 on 2010/12/03 (except for CoRoT-14 which 
was obtained on 2011/01/11) and interpreted using the N2 data 
description document 11 !. The C0R0T data are of two types: short- 
cadence and long-cadence. The total integration times for the two 
cadences are 32 s and 512 s, respectively. Many of the C0R0T TEPs 
have both types of data, and these are treated separately in each 
case. 

The Kepler data used are the public data obtained from the 
Multimission Archive at STScI (MAST0). The Kepler quarter 
(Q0) and quarter 1 (Ql) data were used in the original version of 
this work, but the analysis was revised to include quarter 2 (Q2) 
data when these became available. These data were downloaded 
from MAST on 2011/04/11. As with C0R0T, the Kepler data 
come in short-cadence and long-cadence flavours . The total integra- 
tion times are 58.84 876 s JGilliland et aHl20irj) and 29.4244 min 
(Jenkins et al. 2010b), respectively. 

In order to account for slow variations in the mean flux level 
of the systems in the C0R0T and Kepler data (both astrophysi- 
cal and instrumental), the data surrounding each transit were ex- 
tracted from the full time-series and normalised to unit flux using a 
straight-line fit to the out-of-transit datapoints. In general the num- 
ber of points used was roughly double the number in the transit, 
split equally before and after the transit. A straight line was used 
because it was adequate for the job and it cannot distort transit pro- 
files. In a few cases entire transits were rejected because of a jump 
in observed flux during the transit or because the observations did 
not cover all of the transit event. 

In cases where the number of datapoint was too large for 
Monte Carlo simulations to be completed in a reasonable length 
of time, the datasets were reduced down to a manageable size. This 
was done by converting the timestamps to orbital phase, sorting 
them, and then binning N consecutive datapoints where N ranged 
from 20 to 100. The value of N was carefully chosen to ensure 
the sampling rate was still sufficiently high that numerical integra- 
tion was not needed. The averaging process tends to convert cor- 



8 A few TEPs undergo partial ecl ipses and therefore are technically not 
"transiting" planets (e.g. WASP-34; [Smallev et alj201ll) . 

9 It is technically possible to have two transits per orbit due to gravitational 
lensing in binary systems containing degenerate objects (Marsh 2 0011) but 
no such system is known. 



10 |http : //idoc-corot . ias .u-psud. f r/| 

11 http: //idoc-corotn2 -public .ias.u-psud.fr/jsp/doc/ 
DescriptionN2vl . 3 .pdf 

12 http : //archive . stsci . edu/kepler/ data-search /search .php 
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Table 4. Parameters from the light curve analyses presented here and in previous works, and used here to determine the physical properties of the TEPs. The 
orbital periods are either from this work or from the literature, and the bracketed numbers represent the uncertainty in the preceding digits. Systems for which 
orbital eccentricity or third light was accounted for are indicated with a * in the column marked "e?" or "L3I", respectively. 



System 



Orbital period 
(days) 



L3? Orbital inclination, 
i (degrees) 



Fractional stellar 
radius, r\ 



Fractional planetary 
radius, rj, 



Reference 



C0R0T-I 


1.5089686 (6) 


CoRoT-2 


1.7429935 (10) * 


CoRoT-3 


4.2567994 (35) 


CoRoT-4 


9.20205 (37) 


CoRoT-5 


4.0378962(19) * 


C0R0T-6 


8.886593 (4) 


CoRoT-7 


0.853585 (24) 


C0R0T-8 


6.21229 (3) 


CoRoT-9 


95.2738 (14) * 


CoRoT-10 


13.2406 (2) * 


CoRoT-11 


2.994330(11) 


CoRoT-12 


2.828042 (13) 


CoRoT-13 


4.035190(30) 


CoRoT-14 


1.51214(13) 


CoRoT-15 


3.06036 (3) 


HAT-P-4 


3.0565195(25) 


HAT-P-7 


2.2047304 (24) 


HAT-P-11 


4.88781501 (68) * 


HD 17156 


21.216398(16) * 


HD 80606 


111.4367(4) * 


Kepler-4 


3.213658 (38) 


Kepler-5 


3.548469 (15) 


Kepler-6 


3.2347020 (33) 


Kepler-7 


4.8854948 (82) 


Kepler-8 


3.5225047 (76) 


KOI-428 


6.87349 (64) 


LHS6343 


12.71382 (4) * 


TrES-2 


2.47061323 (7) 


TrES-3 


1.30618700(72) 


WASP-3 


1.8468373 (14) 


WASP-7 


4.9546416 (35) 


XO-4 


4.1250828(40) 



84.42 ±0.31 


0.2073 ± 0.0020 


0.02924 ±0.00041 


87.45 ±0.34 


0.1478 ±0.0023 


0.02462 ± 0.00035 


85.80 ±0.77 


0.1266 ±0.0067 


0.00857 ± 0.00056 


89 96+ - 04 


0.0585±° ;°°« 


n nnfins+o 00060 
U.0UbU8_ 00020 


86.24 ±0.53 


0.0977 ± 0.0067 


0.01 129 ±0.00091 


88.88 ±0.25 


0.0567 ±0.0013 


0.00662 ±0.00019 


79.6 ±3.2 


0.264 ±0.039 


0.0047 ±0.0012 


87.44 ±0.56 


0.0659 ±0.0056 


0.00537 ± 0.00060 


89.97 ±0.13 


0.01083 ±0.00072 


0.001233 ±0.000089 


88.57±0.18 


0.0326 ± 0.0023 


0.00424 ± 0.00038 


83.13±0.19 


0.1452 ±0.0022 


0.01549 ±0.00023 


85.79 ±0.43 


0.1235 ±0.0035 


0.01638 ±0.00074 


85.27 ±0.47 


0.1161 ±0.0053 


0.01 173 ±0.00063 


79.7 ±1.4 


0.206 ±0.019 


0.0181 ±0.0013 


89.9i2.; 

89 - 2 -1.5 


n iqc+0.034 
u.±oo_ 009 

1666+ 00080 
U.lOOO_ 0027 


OOIOQ+O' 0034 
u.urua_ 0008 

01431+0000 72 

U.Ui*ol_ .00028 


83.40 ±0.12 


0.23880 ±0.00095 


0.018401 ±0.000090 


89.36 ±0.36 


0.06148 ±0.00082 


0.003604 ±0.000071 


86.94 ±0.34 


0.04222 ± 0.00079 


0.003107 ±0.000081 


89.232 ±0.029 


0.01056 ±0.00024 


0.001045 ±0.000019 


89.2+.0;S 
87.li ;° 
89.9lJ!; 7 


n 1KQ+0.026 

u.iocs_ 010 

1445+0- 0042 
U.lfiO_ 0033 

1321+ 00017 
u.±ozi_ 0006 


00391+0000 76 

u.uucsai_ 00034 

01163+ 01 :i -' 

U.UllOJ_ 00027 

01259+ 000022 

U.U1Z09_ 00007 


85.31 ±0.43 


0.1494 ±0.0035 


0.01250 ±0.00041 


84.23 ±0.16 


0.1432±0.0018 


0.01360 ±0.00021 


86.5+l 5 6 


130+0.029 

U.lO3_ Q19 


0079+ 00021 

U.UUI»_ 0013 


89.247 ±0.088 


0.02388 ± 0.00090 


0.00489 ±0.00018 


83.925 ±0.030 


0.12568 ±0.00041 


0.015979 ±0.000027 


81.93±0.13 


0.1682 ±0.0014 


0.02750 ±0.00035 


83.72 ±0.39 


0.1994 ±0.0032 


0.02125 ±0.00041 


87.03 ±0.93 


0.1102±0.0061 


0.01053 ±0.00070 


89.9+.°] 


1300+ 00283 
U.lOUU_ 0051 


01124+ 000334 
U.U1±Z4_ 00054 



This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

Southworth et al. (2011) 

Paper III 



related noise into white noise, but does not cause a problem for 
the error estimates of the resulting parameters. The phased and 
binned datapoints will be referred to as 'normal points' below, a 
term which was once used extensively in studies of eclipsing bina- 
ries. At the request of the referee I have analysed the C0R0T 32s- 
sampled data for C0R0T- 1 both with and without phase-binning, in 
order to check that the binning process does not affect the solutions. 
The resulting light curve parameters differ by negligible amounts, 
supporting the use of phase-binning below. 



6 RESULTS FOR INDIVIDUAL SYSTEMS 

In this section I present the photometric (JKTEBOP) and absolute- 
dimensions (JKTABSDIM) analyses of 32 TEPs based on published 
high-quality data obtained from space missions, complemented 
with ground-based data where possible. The results are obtained us- 
ing the same methods as those in Paper III, leading to homogeneous 
measurements for a sample of 58 TEPs. The final JKTEBOP results 
of the 32 TEPs are collected in Table[4] which also includes the 
orbital periods and indicates for which systems a non-circular or- 
bit was adopted. The adopted spectroscopic parameters (T c g , [^1 
and Ka) are given in Table[5] with a lower limit on the errorbars 



of ±50 K for T e a and ±0.05 dex for [§] (see Paper II). Extensive 
tables of results, plus a comparison with literature values, can be 
found in the online-only Appendix. 



6.1 C0R0T-I 

C0R0T-I was the first TEP discovered by the C0R0T satellite 
(Barge et al. 2008), and was originally called CoRoT-Exo- 1. The 
C0R0 T data have been subj ected to TTV anal yses by|B ean (2009) 
and ICsizmadia et alj J2010I) with null results. iGillon et all (2009) 
observed a transit in the i?-band using VLT/FORS2 and an oc- 
cultation at 2.09 ^im using VLT/HAWKI. The transit data are of 
very high precision (0.52 m mag scatter) but do not fully sample 
the transit. Pont et alj ( 120 id) presented the same i?-band data sup- 
plemented with a smaller number of B-band o bservations of the 
same transit taken with the same instrument. IPont et alj q2010h 
also obtained RV observations during a transit and detected a 
Rossiter-McLaughlin (RM) effect which show that the orbital axis 
of CoRoT-lb is not aligned with t he st ellar spin axis. Occ upa- 
tion studies bylAlonso et alj J2009bl) and lRogers et alj J2009h and 
iDeming et alj ( 120 111) have found no evidence for orbital eccentric- 
ity, so I assumed that the orbit is circular. 

The C0R0T data show 20 transits observed at long cadence 
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Table 5. Measured quantities for the parent stars which were adopted in the analysis presented in this work. 



System 


Velocity amplitude ( m s 1 ) 


Toff (K) 


Reference 




"Fc" 
_ H . 


Reference 


CoRoT-1 


188 ±11 


Barge et al. (2008) 


5950 ±150 


Barge et al. (2008) 


-0.30 ±0.25 


Barge et al. (2008) 


CoRoT-2 


603 ±18 


Gillon et al. (2010a) 


5696 ± 70 


Chavero et al. (2010) 


0.03 ±0.06 


Chavero et al. (2010) 


CoRoT-3 


2170 ±27 


Triaud et al. (2009) 


6740 ±140 


Deleuil et al. (2008) 


-0.02 ±0.06 


Deleuil et al. (2008) 


CoRoT-4 


63 ±6 


Aigrain et al. (2008) 


6190 ±60 


Moutou et al. (2008) 


0.05 ±0.07 


Moutou et al. (2008) 


CoRoT-5 


59.1+J;? 


Rauer et al. (2009) 


6100 ±95 


Rauer et al. (2009) 


-0.25 ±0.06 


Rauer et al. (2009) 


CoRoT-6 


280 ± 30 

5.04 ±1.09 


Fridlund et al. (2010) 


6090 ± 70 
5250 ± 60 


Fridlund etal. (2010) 


-0.20 ±0.10 
0.12 ±0.06 


Fridlund etal. (2010) 


CoRoT-7 


Hatzes et al. (2010) 


Bruntt et al. (2010) 


Bruntt et al. (2010) 


CoRoT-8 


26±4 


Borde et al. (2010) 


5080 ± 80 


Borde etal. (2010) 


0.31 ±0.05 


Borde etal. (2010) 


CoRoT-9 


38 ±3 


Deeg et al. (2010) 


5625 ± 80 


Deeg etal. (2010) 


-0.01 ±0.06 


Deeg etal. (2010) 


CoRoT-10 


301 ± 10 


Bonomo et al. (2010) 


5075 ± 75 


Bonomo et al. (2010) 


0.26 ±0.07 


Bonomo et al. (2010) 


CoRoT-11 


280 ±40 


Gandolfi et al. (2010) 


6440 ±120 


Gandolfi et al. (2010) 


-0.03 ±0.08 


Gandolfi et al. (2010) 


CoRoT-12 


125.5+*-° 


Gillon etal. (2010b) 


5675 ± 80 


Gillon et al. (2010b) 


0.16±0.10 


Gillon et al. (2010b) 


CoRoT-13 


157.8 ±7.7 
1230 ± 34 


Cabrera etal. (2010) 
Tinglevetal. (2011) 


5945 ± 90 
6035 ± 100 


Cabrera etal. (2010) 
Tinglevetal. (2011) 


0.01 ±0.07 
0.05 ±0.15 


Cabrera etal. (2010) 


CoRoT-14 


Tinglevetal. (2011) 


CoRoT-15 


7360 ±110 


Bouchvetal. (2011) 


6350 ± 200 


Bouchvetal. (2011) 


0.1 ±0.2 


Bouchvetal. (2011) 


HAT-P-4 


81.1 ±1.9 


Kovacs et al. (2007) 


5860 ± 80 


Kovacs et al. (2007) 


0.24 ±0.08 


Kovacs et al. (2007) 


HAT-P-7 


211.8±2.6 


Winn et al. (2009a) 


6350 ± 80 


Pal et al. (2008) 


0.26 ±0.08 


Pal et al. (2008) 


HAT-P-11 


11.8±0.9 


Hiranoetal. (2011) 


4780 ± 50 


Bakos et al. (2010) 


0.31 ±0.05 


Bakos et al. (2010) 


HD 17156 


272.7 ±2.1 


Winn et al. (2009c) 


6079 ± 56 


Fischer et al. (2007) 


0.24 ±0.05 


Fischer et al. (2007) 


HD 80606 


476.1 ±2.2 


Winn et al. (2009b) 


5574 ± 72 


Santos et al. (2004) 


0.34 ±0.05 


Gonzalez et al. (2010) 


Kepler-4 


**i\\ 


Borucki etal. (2010a) 


5857 ± 120 


Borucki et al. (2010a) 


0.17 ±0.06 


Borucki et al. (2010a) 


Kepler-5 


227.5 ±2.8 


Koch etal. (2010) 


6297 ± 60 


Koch et al. (2010) 


0.04 ±0.06 


Koch et al. (2010) 


Kepler-6 


80.9 ±2.6 


Dunham et al. (2010) 


5647 ± 50 


Dunham etal. (2010) 


0.34 ±0.05 


Dunham etal. (2010) 


Kepler-7 


42.9 ±3.5 


Latham et al. (2010) 


5933 ± 50 


Latham etal. (2010) 


0.11 ±0.05 


Latham etal. (2010) 


Kepler-8 


68.4 ±12.0 


Jenkins et al. (2010a) 


6213 ±150 


Jenkins et al. (2010a) 


-0.055 ±0.05 


Jenkins et al. (2010a) 


KOI-428 


179 ±27 


Santerne et al. (2011) 


6510 ±100 


Santerne etal. (2011) 


10+ 015 
u - ±u -o.io 


Santerne et al. (2011) 


LHS 6343 


9600 ± 300 


Johnson etal. (2011) 


3300 ± 200 


This work 


0.04 ± 0.08 


Johnson etal. (2011) 


TrES-2 


181.3±2.6 


O'Donovan et al. (2006) 


5850 ±50 


Sozzetti et al. (2007) 


-0.15±0.10 


Sozzetti et al. (2007) 


TrES-3 


369 ± 11 


Sozzetti et al. (2009) 


5650 ± 75 


Sozzetti et al. (2009) 


-0.19 ±0.08 


Sozzetti et al. (2009) 


WASP-3 


286.5 ±7.8 
97 ±13 


This work 

Hellier et al. (2009b) 


6400 ± 100 
6400 ±100 


Pollacco et al. (2008) 


0.00 ± 0.20 
0.00 ±0.10 


Pollacco et al. (2008) 


WASP-7 


Hellier et al. (2009b) 


Hellier et al. (2009b) 


XO-4 


165.8 ±6.2 


Narita etal. (2010) 


6397 ± 70 


McCullough et al. (2008) 


-0.04 ±0.05 


McCullough et al. (2008) 



(512 s) and 17 observed at short cadence (32 s). Each transit was 
normalised (Sect.[5]l and the two cadences were treated separately 
in the JKTEBOP analysis. After a preliminary fit the 32 s data were 
phase-binned with each phased datapoint representing 50 original 
ones. Seven out of 546 of the 512 s datapoints were rejected as 4<r 
outliers. The 512 s data were modelled using numerical integration 
with ATjnt = 3. The i?-band and B-band data were also studied 
(FigJJ. 

The results for the CoRoT 32 s data are given in Table A 1 
and are of sufficiently high quality for the LD-fitted solutions 
to be adopted. LD-fit/fix solutions were adopted for the CoRoT 
512 s (Table A2) and the VLT i?-band (Table A3) data, and in both 
cases correlated noise was found to be unimportant. Correlated 
noise is moderately important for the B-band data (the residual- 
permutation errorbars are larger than the Monte-Carlo errorbars) 
but the data were good enough to allow the LD-fit/fix solutions to 
be adopted (Table A4). 

The photometric results are given in Table A5 and show good 
agreement except for k (which has a reduced \ 2 of xS — 3.2). This 
phenomenon has been n oted several times before (Paper I, Paper II, 
ISouthworth et alj|2009bl) and is attributable to correlated residuals 
(both instrumental and astrophysical) in the data. The B-band data 
cannot match their counterparts so I combine the other three sets of 
results together to get the final photometric parameter s. These are 
in go od agreement with those of the discovery paper (Barge et al. 
l2008h and of lBeanl J2009t) . but not with lGillon et al] te009l) 



The physical properties of CoRoT- 1 have been calculated us- 
ing JKTABSDIM. The individual solutions are given in Table A6 
and compared with litera ture values, where a good agreement is 
found (except for px from Gillon et al. 2009 due to their somewhat 
different photometric results). My measurements of <7b, which is 
quite similar to the Earth's surface gravity, and Q appear to be the 
first ones published. Our understanding of the system would benefit 
from further spectroscopy to give improved values of T B a , Ka and 
[??] in particular. 

6.2 CoRoT-2 



The discovery of CoRoT-2 was announced bv lAlonso et alj (2008) 
based on a 142 d light curve from the CoRoT satellite which ex- 
hibited transit events every 1.74 d and significant starspot activity 
on longer timescales. CoRoT-2A is a young st a r hos ting a rel- 
atively massive planet (3.6M.t„ p ). iBouchv et al] 120081) obtained 
RV measurements through transit and modelled the RM effect to 
find that the planetary orbital and projected stellar spin axes are 
aligned to within 7.2° (1 .6a ). The occul tation of the planet h as been 
found in the CoRoT data bv lAlonso et a l. ( 2009a) and S nellen et al] 
d2010b|). and its time of o ccurrence is consistent with a circular 
orbit. iGillon et al] ( I2010al) presented Spitzer observations of the 
secondary eclipse at 4.5 pm. and 8 pm, finding a small but signifi- 
cant eccentricity characterised as e cos 10 — — 0.00291+oQQQgi and 
esinoj = 0.0139l~ 0OOg4 . Similar results have been obtained by 
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Figure 5. Phased light curves of CoRoT- 1 compared to the best fit found 
using JKTEBOP and the quadratic LD law. The residuals are plotted at the 
base of the figure, offset from unity. The purple line through some data 
show the best-fitting model without numerical integration - in these cases 
the difference between this model with and without numerical integration 
is shown by another purple line through the residuals. From top to bottom 
the light curves are the binned CoRoT 32 s data, the CoRoT 512 s data, the 
FORS2 i?-band data and the FORS2 B-band data. 



ity darkening. In the case of starspots this is certainly not a reli- 
able assumption. However, the effect will average out over a large 
number of transits if the starspots do not show a preference for 
particular latitudes. But we know from the Sun that this is not 
the case: sunspots appear mostly within 30° of the equator, and 
their preferred latitudes vary throughout the 1 1 yr solar activity cy- 



cla I A different effect is noticed in the study of active stars (e.g. 
lOlah et alJI 19971 : [Barnes 2005) and eclipsing bi naries, where li ght 
curve solutions often favour large polar spots (Hilditch 200 lh. It 
is not usually possible to account for this problem in TEP stud- 
ies because of the difficulty in detecting starspots outside the area 
eclipsed by the planet. The effect of star spots on the analys is of 
the CoRoT-2 syst em has been studied bv lCzeslaet al. (2009) and 
iHuber et al.l J201fj|) , who agree that the spot coverage on the chord 
of the planet transit is greater than the average for the stellar disc, 
and that the planet is therefore a few percent larger than standard 
analyses would suggest. 



In the current work I have modelled the CoRoT 32 s light 
curve, which covers 79 transits. The 512 s data are not used as they 
only spread over three transits with partial coverage. 148 points 
of the 32 s data were rejected by a Aa clip and the remaining 
50 666 points were phase-binned by a factor of 200 into 254 nor- 
mal points. In the JKTEBOP analysis I did not account for starspots, 
as this is beyond the scope of the current work, so have produced 
a baseline solution which is more easily comparable to the results 
for other TEPs. This is equivalent to assuming that the starspots 
affect all parts of t he star equally on average. A third light of 
L 3 = 0.053±0.003 JAlonso et alj2008h and the small eccentricity 
JGillon et alJuOlOal) is acc ounted for in the wa ys described in Pa- 
per III. The VLT data from lGillon et alj < !2010ah were not modelled 
because they only cover half a transit. 
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Figure 6. The 32 s light curve of CoRoT-2 compared to the JKTEBOP best 
fit. Other comments are the same as Fig. [5] 



iDeming et al.l J2011I) from Warm-Spitzer observations at 3.6 /im: 
ecosw = -0.0030 ± 0.0004. 

The spot activity of CoRoT-2 A deserves mention. It is an ac- 
tive star with a rotation period of only 4.5 d and probable differ- 
ential rotation JLanza et alj[2009: Froh lich etai¥2 009: Hube r et al] 



l201fJ : ISilva-Valio & L anza 2011). One of the standard assumptions 
of modelling transit light curves of TEPs is that the surface bright- 
ness of those parts of the star eclipsed by the planet is the same 
as that of the rest of the star, modulo effects such as LD and grav- 



For the phase-binned 32 s data I found that correlated noise is 
not important and that the LD-fitted solutions are reliable enough to 
be adopted as the final photometric results (Table A7). A compar- 
ison with literature measurements (Table A8) shows an acceptable 
agreement except compared to the more sophisticated analysis of 
ICzesla et ail J2009l) . My best fit to the 32 s data is plotted in Fig.[6] 



The young age of CoRoT-2 (it is expected to be of Pleiades 
age) manifests itself in the JKTABSDIM analysis by the best solu- 
tions in some cases being on the zero-age edge of the grids of theo- 
retical predictions. A modestly lower T s g would alleviate this prob- 
lem whilst not having a significant effect on the resulting physical 
properties. The edge effects can lead to underestimated errorbars 
for the affected solutions (see Paper II) but this is not a problem 
because the final results then rest on the errorbars from the unaf- 
fected solutions. The edge effect do, though, cause an increase in 
the systematic errorbars (Table A9). 



CoRoT-2A would benefit from an improved T c s measure- 
ment: the discovery value of 5625 ± 120 K JAlonso et alj|2008h 
is comp aratively uncertain, the alterna tive measurement of 5608 ± 
37 K bv lAmmler-von Eiff et al.l J2009I) accompanies a log g which 

is too high (4.71 versus 4.53), and the value of 5696 ± 70 K found 

I — 1 J — •!, 

bv lChavero et al. (2010) is higher so causes stronger edge effects. 

A more precise measurement of K\ would be useful. 
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Figure 7. The 32 s (upper) and 512 s (lower) CoRoT light curves of CoRoT- 
3. Other comments are the same as Fig. [5] 
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Figure 8. The 32 s (upper) and 512s (lower) CoRoT light curves of CoRoT- 
4. Other comments are the same as Fig.|5] 



6.3 CoRoT-3 

CoRoT-3 was announced bv lDeleuil et alj ( 120081 ) and was the first 
brown-dwarf-mass object with a precisely measured mass and ra- 
dius. Its F3 V host star is rather hot (6740 K), making photometry 
and RV measurements difficult due to a shallow transit and high 
stellar v sin i. An RM study has been presented by Triaud et al. 
(2009), whose measurements are consistent with a circular and ax- 
ially aligned orbit. 

The 32 s data cover 19 transits with a total of 196 691 data- 
points. Of the 19046 points in the regions of the transits, 20 were 
rejected by a 4a clip and the remainder were phase-binned by a 
factor of 40 to get 477 normal points. The 512 s data cover 16 tran- 
sits and were fitted using 7Vi nt = 3. A circular orbit was assumed 
and a third light of L 3 = 0.082 ± 0.007 was adopted JDeleuil et all 
l2008h . 

The results for the 32 s data are given in Table A10: corre- 
lated noise is not important (as usual for phase-binned datapoints) 
and the LD-fitted solutions yield the lowest scatter and reason- 
able LDCs. For the 512 s data (Table All) correlated noise is again 
unimportant and the LD-fit/fix solutions are best. The best fits are 
plotted in Fig.|7] The final photometric parameters are the weighted 
means of the 32 s and 512 s values. They are compared to published 
values in Table A 12, where a very good agreement is found. 

The measured physical properties of CoRoT-3 are given in Ta- 
ble A13 and agree well with literature values. I provide the first 
measurement of T c ' q and O. The star could do with a better T B g 
measurement, and an improved light curve would also be useful. 



6.4 CoRoT-4 



CoRo T-4 was discovered bv lAigrain et al. ( 2008) and lMoutou et al.1 
(2008). It has a relatively long orbital period of 9.2 d, and thus 
one of the lower T cq s despite having a late-F host star. There are 
three transits each in the 32 s and the 512 s data. The 32 s data were 



13 The spot characteristics of the Sun are nicely captured in the 
Maunder butterfly diagram. An up-to-date version can be found at: 

|http: //s olarscience . msf c .nasa . gov / image s/bfl y .gif| 



phase-binned by a factor of 10 after 4a clipping to remove a small 
number of outliers. The 512 s data were modelled using JVj n t = 3. 

The JKTEBOP solutions favour a central transit, which causes 
the photometric parameters to have asymmetric errorbars. Corre- 
lated noise is unimportant for the phase-binned 32 s data but is im- 
portant for the 512 s data. In both cases the LD-fit/fix results are the 
best (Tables A14 and A15). These were combined into final pho- 
tometric parameters by multiplying their probability density func- 
tions (Table A 16). The best fits a re shown in F ig.[H The agreement 
between my results and those of I Ai grain et alj ( 120081) is excellent. 

The physical properties of the CoRoT-4 syste m are given in 
Table A17 and show good agreement with those of iMoutou et al.1 
(2008) except for larger errorbars for some properties (notably Ma 
and Ra)- I provide the first published measurements of pa, gb and 
0. The long orbital period and short observing run conspire to- 
gether to allow only six transits to be observed by CoRoT: a better 
light curve would be useful, as would a more precise Ka- 



6.5 CoRoT-5 

Discovered bv lRauer et al.1 (J2009), CoRoT-5 is a fairly normal sys- 
tem containing a low-density TER The 32 s data contain 23 tran- 
sits, which were cut out of the light curve, 4a clipped and phase- 
binned as usual. The 512 s light curve contains six transits and these 
data were solved using iVj n t = 3. iRauer et a l. (2009) found a pre- 
liminary third light of 8.6% (no uncertainty) so I ad opted L3 = 
0.086 ± 0.020 for the JKTEBOP analysis. Note that IRauer et al.1 
(2009) removed this contaminating light from the light curve prior 
to modelling it, so neglected the uncertainty in L3. Eccentricity is 
significant at the 3a level so this was included via the constraints 



-0.057t ; 2 andesinw = 



-0.071 



+0.147 



The JKTEBOP results are shown in Tables A18 and A19, and 
the best fits are plotted in Fig.|9] In both cases correlated noise was 
unimportant and the LD-fit/fix solutions were adopted. The two 
sets of results do not appear to agree very well, but are consistent 
to within la. The errorbars are large primarily due to the poorly 
constrained orbital parameters: ecosu and esino; are both uncer- 
tain and the latter is significantly correlated with ta and n>. The 
combined results (Table A20) are in poor agreement with those of 
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Figure 9. The 32 s (upper) and 512 s (lower) CoRoT light curves of CoRoT- 
5. Other comments are the same as Fig. [5] 
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Figure 10. The 32 s light curve of CoRoT-6. See Fig.|5]for details. 



iRauer et al.l J2009T) . presumably due to differences in the analysis 
methods. 

The JKTABSDIM results are given in Table A21 and show a 
good agreement between the different model sets (and with the 
dEB constraint discussed in Sect. l3.2t . As expected given the dis- 
crepant light curve results, I fin d physical properties which differ 
from those of lRauer et all ( 120091) by up to 2a. I find star and plane t 
radii which are notably smaller than those of lRauer et alj (2009), 
and errorbars which are substantially larger (by a factor of 5 for 
Ma). The properties of CoRoT-5 are more uncertain than previ- 
ously thought, and additional RV measurements are the best route 
to fixing this. 



6.6 CoRoT-6 

Comparatively speaking, CoRoT-6 is a massive (2.96 Mj up ) and 
dense (1.66 pj up ) TEP with a long orbital period (8.89 d) around a 
moderately active F9 V star. The discovery light curve from CoRoT 
JFridlund et al.l2010h contains 33 1 397 datapoints at short cadence, 
covering 14 c omplete transits. A study of the starspot activity has 
been given bv lLanza et alj ( 120111) . 
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Figure 11. The 32 s light curve of CoRoT-7. See Fig.|5]for details. 



I adopted e = and L- A = 0.028 ± 0.007 JFridlund et ail 

I2OIOI) . 40 out of the 23 802 near-transit datapoints were rejected 
by a 4<r clip and the remainder were phase-binned by a factor of 
50 into 476 normal points. The JKTEBOP results a re given in Ta- 
ble A2 2 and are in good agreement with those of iFridlund et alj 
OO 10). Correlated noise is completely negligible and the LD-fit/fix 
solution is the best choice (Table A23). The best fit is plotted in 
Fig.Q! 

The physical properties of C0R0T-6 are given in Table A24, 
where the agree ment between theore tical models, the dEB con- 
straint and with IFridlund et alj yOlOJ) is excellent. I provide the 
first measurements of pA, <?b and O. The system could do with 
more spectroscopy to obtain an improved Ka and [|^1 . 

6.7 CoRoT-7 

CoRoT-7 is one of the most important known planets. At the time 
of discovery it was both the smallest and least massive TEP discov- 
ered so far. The CoRoT light curve shows a transit of depth only 
0.034% recurring on a period of only 0.85 d. A detailed analysis of 
these data, plus t wo high-precision RVs, proved the planetary na- 
ture of this object dLeger et al. 2009). Extensive RV measurements 
using the HARPS spectrograph dOueloz et alj|2009l) confirmed the 
nature of the planet, with amass of 4.8 ± 1.8 M$, and allowed the 
discovery of a second one which is more massive (8.4 M$) and on 
a longer-period orb it (3.70 d). A study of the atmospheric param- 
eters of the star by iBruntt et alj 020101) yielded a smaller Ra and 
therefore Rh- 

CoRoT-7 A is a G9 dwarf star with significant chromospheric 
activity. This activity causes larger RV variations that those of 
the planets' o rbits, making the pr operties of the system somewhat 
controversial. lLanza et al.l (120100 studied the stellar activity and 
starspots and found a false-alarm probability of <10~ 4 that the 
RV oscillations attributed to CoR oT-7b and CoRoT-7 c are spu- 
rious effects of noise and activity. lHatzes et alj 020100 performed 
a comprehensive re-analysis of the RVs, finding a larger mass for 
CoRoT-7b (6.9 ± 1.4M©), confirming the RV signal of the sec- 
ond planet (CoRoT-7c) and tentatively detecting a third RV sig- 
nal which could be caused by another planet, CoRoT- 7d, with a 
mass of 16.7 ±0.4 Me and a period of 9.02 d. However. IPont et alj 
d201ll) performed a similar analysis, including a phenomenological 
model to describe the properties and evolution of many starspots, 
and found the RV signal of CoRoT-7b to be a lot smaller and 
only significant at the 1.2a level (M h = 2.3 ± 1.5M©). The 
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Figure 12. The 32 s (upper) and 512 s (lower) CoRoT light curves of 
CoRoT-8. Other comments are the same as Fig.[5] 



huge number of parameters in their starspot model could be ex- 
pected to lead to more hazy r e sults c ompared to those of other re- 
searchers. iFerraz-Mello et al.l J2010I) showed that previous analy- 
ses have tended to remove some of the planetary RV signal when 
squashing the effects of the stellar activity, and thus underestimate 
the masses. They found Mb = 8.5 ± 1.5 Me and M c sin i c = 
13.5 ± 1.5 M®, an d that CoRoT-7 d is an artefact rather than an as- 
trophysical signal. iBoisse et alj ( 120111) also modelled the starspot- 
induced RV variations and confirmed planets b and c with false- 
alarm probabilities of < 5 x 10~ 4 . They found masses of Mb = 

13.2 ± 4.1 M ffi . Just before the 



5.7 ± 2.5 M and M c shr 5 i c 
present paper was submitted, lHatzes et al.1 J20 1 lh produced a re- 
analysis which confirmed their previous results for CoRoT-7 b. The 
difference in Ka (5.15 ± 0.95 m s _1 versus 5.06 ± 1.06ms -1 ) is 
so small that there was no need to recalculate my results. 

The CoRoT light curve contains 308 947 datapoints, all at 
32 s cadence. The 47 702 points near transits were phase-binned 
by a factor of 300 to give 160 normal points. The transit is so 
shallow (Fig.lllb that there was no point in calculating LD-fitted 
models. The results for the LD-fit/fix solutions (Table A25) show 
a good mutual agreement whereas the LD-fixed soluti ons do not. I 
theref ore adopted the LD-fit/fix solutions. Compared to lLeger et al.l 
(2009), my results are in good agreement but have much larger er- 
rorbars (Table A26). This wraps over into the physical properties 
(Table A27). My analysis approach may not be optimal for this sys- 
tem, as I do not apply any external constraints to the light curve 
model (e.g. and spectroscopically derived log g or pa), but chang- 
ing this would cause an inhomogeneity with the other TEPs treated 
in this work. CoRoT-7 would benefit from improved photometric 
and RV observations, although these are observationally highly de- 
manding. 



6.8 CoRoT-8 

CoRoT-8 was discovered bv lBorde et al.l a201Cl) and is a small and 
low-mass planet in a relatively long-pe riod orbit (6.21 d) around a 
metal-rich K dwarf. lBorde et alj <20101) quote a third light value of 
0.9% without an errorbar; I adopted L 3 — 0.009 ± 0.003 in my 
photometric analysis. The 32 s data cover 12 transits with 19413 



datapoints. 46 were rejected by a 4a clip and the remainder were 
phase-binned by a factor of 40 to give 483 normal points. The 
512 s observations harbour 1 1 transits; nine out of 946 points were 
rejected by a 4<r clip and the remaining data were solved using 
Nut = 3. 

In the course of extracting the tr ansits from the Co RoT data 
it became clear that the ephemeris in iBorde et alj (201fj) predicts 
transits to occur too early. I therefore binned up the 32 s data by a 
factor of 16 to match the 512 s data and fitted both datasets together 
to get a new orbital ephemeris: 

To = HJD 2 454 239.03311(78) + 6.212381(57) x E 

where E is the number of orbital cycles after the reference epoch 
and the bracketed quantities denote the uncertainty in the final digit 
of the preceding number. The errobars come from Monte Carlo and 
residual permutation simulation s, which are in g ood agreement. 
Compared to the ephemeris of IBorde et alj < 20 101) I obtain larger 
errorbars, a consistent orbital period, and a To which is later by 
0.059 d (67<j). P. Borde (private communication) has kindly con- 
firmed that this discrepancy has been noted elsewhere, by two ama- 
teur astronomers who have found that the transits of CoRoT-8 occur 
later than predicted from the ephemeris in the discovery paper. 

The 32 s and 512 s datasets were then fitted individually, the 
latter with iVint = 3. For both, correlated noise is unimportant 
and the LD-fit/fix solutions are adopted (Tables A28 and A29). 
The best fits are shown in Fig.[T2j The combined solution of the 
two dataset s (Tab le A30) does not agree well with the results of 
IBorde et al.l 12010J), in particular t a (1.6<t) and k (2.8<r). The error- 
bars found bv lBorde et alj (2010) seem to be too small (e.g. 0.1° 
compared to my 0.6° for i), especially in light of the disagreement. 

The physical properties of CoRoT-8 were not straightforward 
to derive. The long evolutionary timescale of the host star caused 
the best fits for the Teramo and VRSS models to be slightly dis- 
crepant with ages formally in excess of a Hubble time. How- 
ever, reasonable solutions could be found by restricting the age 
ra nge to 0-10 Gyr. My results (Table A3 1) again differ from those 
of IBorde et alj J2010I) , most obviously for the radii of the star 
(0.90 ± 0.09 versus 0.77 ± 0.02) and planet (0.71 ± 0.08 ver- 
sus 0.57 ± 0.02). The masses of the components agree well with 
my results, and I present the first measurements of (?b, T c ' q and O. 
I find that the density of the planet is less than half the value in the 
discovery paper. Further light curves and RVs would be useful for 
this object, in order to improve the system parameters and confirm 
the orbital ephemeris determined above. 



6.9 CoRoT-9 

Th e discovery ofthe 93.7-d period CoRoT-9 system was announced 
bv lDeeg et a l. (2010): the system has the longest P or b of all TEPs 
bar HD 80606 (which was discovered in the course of an RV rather 
than a photometric survey). The CoRoT data cover only two tran- 
sits of which one lacks coverage of ingress. The data were cut 
from the full light curve and phase-binned by a factor of 5 to yield 
946 datapoints (Fig.ll3t. The orbita l shape was constrained using 
e = 0.11 ± 0.04 and u = 37±! 7 foeeg et alj|20ich . The JKTE- 
BOP solutions show that correlated noise is not important and that 
the LD-fit/fi x alternative is be st (Table A32). The agreement with 
the results of lDeeg et alj J20101 ) is excellent, both for the photomet- 
ric parameters (Table A33) and the physical properties (Table A34). 
CoRoT-9 would benefit from further light curves and RVs. 
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Figure 13. The 32 s light curve of CoRoT-9. See Fig.[5]for details. 
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Figure 14. The CoRoT 32s (upper) and 512s (lower) light curves of CoRoT- 
10. Other comments are the same as Fig.[3] 



6.10 CoRoT-10 



CoRoT-10 was discovered bv lBonomo et alj J2010I) and comprises 
a relatively massive TEP (2.78 Mj up ) in a long-period (13.24 d) 
and highly eccentric (e = 0.53 ± 0.04) orbit around a Kl dwarf. 
Five transits were observed at 512 s sampling and five more at 32 s 
sampling, of which one was ignored because it is affected by a 
data jump. A third l ight of Ls = 0.055 ± 0.003 was found by 
lBonomoet~aT] ( l2010h . 

The 32 s data have 7190 points in the vicinity of a transit, of 
which 16 were rejected by a 4a clip and the rest phase-binned by a 
factor of 20 to give 369 normal points. No clipping or binning was 
needed for the 512 s data, which were solved using Mnt = 3. The 
best fits are shown in Fig.[l4]and the solutions arranged in Tables 
A35 and A36. M y final photome t ric pa rameters (Table A37) agree 
well with those of Bonomo et alj ( I2010h . 

The physical properties of the system are not precisely defined 
(Table A38). It appears to be a young system, resulting in edge ef- 
fects with the model grids (see CoRoT-2) and therefore larger sys- 
tematic errorbars th an is typical, especiall y for A/a- However, the 
agreement with the iBonomo et alj ( 120 101) parameters is good. A 
better light curve would be beneficial. 
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Figure 15. The CoRoT 32 s (upper) and 512 s (lower) light curves of 
CoRoT-1 1. Other comments are the same as Fig.|3] 



6.11 CoRoT-11 

The star in the CoRoT-1 1 system is one of the earliest- type (F6 V) 
and most rapidly ro tating (40 ± 5kms _1 ) TEP hosts known 
( IGandolfi et alJuOlOO . The planet itself is comparatively massive 
(2.34 Mjup) and has a high T c ' q (1735 K). IGandolfi et alj OoToh 
found a third light of L3 = 0.130 ± 0.015 which is included in the 
JKTEBOP analysis. 

The CoRoT 32 s data have 26 382 points near transits, of 
which 46 were rejected by a 4a clip and the rest phase-binned 
by a factor of 50 to get 527 normal points. The 512 s data con- 
tain 707 near transit of which 4 were rejected by a 4a clip and 
the remainder modelled using 7Vi nt = 3. For both light curves 
the LD-fit/fix solutions are the best and correlated noise is incon- 
sequential (Tables A39 and A40). The fitted data are exhibited in 
Fig-E] The results for the two datasets are unusual in that k is the 
most consistent, whilst the other parameters agree to within less 
good 1.1a. They were combined with this trifling disagreement ac- 
counted for_jn_fhe_ejrrorbars. The final results are very similar to 
those of lGandolfi et alj ( 120100 . and are among the better results for 
the known TEPs (TableA41). This is helped by the relatively low 
i, which means the light curve fits are well-constrained. 

The JKTABSDIM results for different model sets and for the 
dEB constrain t (Tabl e A42) agree well with each other and with 
IGandolfi et alj 1 12010]) . Further spectroscopic study of CoRoT-11 
would be profitable. 



6.12 CoRoT-12 



The discovery paper for the CoRoT-12 system JGillon et alj2010bl) 
presents a system which is rather typical of the known population. 
The 32 s data comprise 242 558 datapoint of which 29 1 14 are ad- 
jacent to one of the 36 transits. A 4a clip discards 76 points and the 
remainder end up in 291 phase-binned normal points. The 512 s 
data cover 11 transits of which 413 points remain after rejecting 
2780 which are away from eclipse and six which are over 4a away 
from a preliminary fit. 

For the JKTEBOP fits I adopted a cir c ular o rbit because the 
e cos uj and e sin u; values in iGillon et alj ( feOlObh are consistent 
with zero at the la level. A third light of L 3 = 0.033 ± 0.005 
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Figure 16. The CoRoT 32 s (upper) and 512 s (lower) light curves of 
CoRoT-12. Other comments are the same as Fig.[5] 
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Figure 17. The CoRoT 32 s (upper) and 512 s (lower) light curves of 
CoRoT-13. Other comments are the same as Fig.|5] 



was used and ./Vint = 3 was employed for the 512 s data. The 
best fits are shown in Fig.[T6]and the model solutions in Tables 
A43 and A44. In both cases correlated noise is unimportant and 
the LD-fit/fix solutions the most reliable. The results (Table A45) 
again show a slight disagreement between the two cadences, with 
k the worst offender (1.7<r) and the other parameters divergent by 
an acceptable 1.2<j. The combine d results agree reasonably well 
with those of iGillon et alj ( I2010bl) . A better light curve would be 
useful. 

The fin a l physic al properties (Table A46) agree very well with 
IGillon et alj J2010b0 . Unusually, many of my error estimates are 
smaller than those found in the literature, which is attributable to 
the standard assumption of e = which neglects any uncertainty 
in the orbital shape. CoRoT-12 could do with further spectroscopic 
observations for both RV and spectral synthesis studies. 



6.13 CoRoT-13 

The discovery paper of this TEP is ICabrera et alj a2010f) . CoRoT 

observed 217 186 datapoints at 32 s cadence, covering 23 transits. 
Of the 22 020 points near transit, 54 were 4<r clipped and the rest 
phase-binned by a factor of 50 to get 441 normal points. Of the 
2403 512 s datapoints, 196 are near transit and were solved using 
Mnt = 3. The quality of the 512 s light curve is relatively poor due 
to the faintness of the star. A third light of L3 = 0.11 ± 0.01 was 
adopted from the discovery paper. 

The best fits of the 32 s and the 512s data are strikingly differ- 
ent (Fig.|17t; I have verified that they arrived in the same datafile. 
The solutions for the 32 s data are in Table A47 and show that cor- 
related noise is unimportant and that LD-fit/fix is to be preferred. 
Models with LDCs fixed at theoretical values have stronger LD and 
return r\ and rb smaller by la. 

The 512 s light curve was difficult to model and suffered from 
instability of solution. This can be seen in most clearly by inspect- 
ing the inclination values in Table A48, which veer from grazing 
in the quadratic LD-fit/fix and square-root LD-fit/fix solutions to 
equatorial in the other 12 solutions. Other parameters are similarly 
affected. A greater consistency might be obtained by guiding the 
offending best fits towards the area of parameter space inhabited 
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Figure 18. Plot of the masses and radii of the known TEPs (dark blue 
crosses for those studied in this series of papers and light red crosses 
for other objects). T he black open circle shows CoRoT-13b (results from 
ICabrera et al. 2010) and the black filled circle shows CoRoT-13b (results 
from this work). The grey dotted lines show the loci where density is equal 
to p.j up , 0.5p Jup and 0.25 pjup. 



by the other best fits, but at the expense of mathematical rigour. A 
better idea is to reject the 512 s data, which are in any case of much 
lower weight than the much more extensive and better-sampled 32 s 
data. For my final photometric parameters I accordingly adopt the 
32 s LD-fit/fix solution. Correlated noise is not important for this 
light curve. 

A comparison with the results of lCabrera et alj 1 120 101) is given 
in Table A49 and shows a poor agreement, most likely because the 
solution in the discovery paper rests on the unreliable 512 s data as 
well as the reliable 32 s data. My results are therefore to be pre- 
ferred, and also have larger and more representative errorbars. 

The final physical properties of CoRoT-13 (Table A50) unsur- 
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Figure 19. The CoRoT 32 s (upper) and 512 s (lower) light curves of 
CoRoT-14. Other comments are the same as Fig.[5] 



prisingly do not agree well with those of ICabreraet alj 12010]) . I 
find radii which are much larger: 26% (3.4<r) for Ra and 41% 
(4.9<t) for _Rb- In turn, g^ and pb are smaller by over a factor of 2. 
The log g\ also decreases from 4.46 ± 0.05 to 4.26 ± 0.04, and the 
spectroscopically derived log g for the host star (4.30 ± 0.10) is in 
slightly better agreement with my value. A comparison of the dif- 
ferent results for CoRoT-13 b in the context of other TEPs is given 
in Fig.Q~8]and shows that my results put it in a regio n of parame- 
ter spa ce occupied by many other TEPs whereas the lCabrera et alj 
12010) results make it unusually dense. An improved light curve 
for the system is urgently needed to verify the revised parameters 
that I find for the system. 



6.14 CoRoT-14 



Discovered bv lTinglev et alj ( 120111) , CoRoT-14 contains a very hot 
Jupiter (T c ' q = 1936 K) on a short-period orbit around an F9 V star. 
The planet's large mass (Mb = 7.7 Mj u p ) is unusual and is second 
only to WASP- 18 faellier et alj|2009al : ISouthworth et al.l J2009ch 
for objects with P or b < 2d. It is rather faint (V — 16.0) so the 
light curve is quite scattered, but the RVs are quite sufficient due to 
its large mass. 

The 32 s data total 217 262 points and 26 279 of these cover 60 
transits; after 4<r-clipping 49 points and phase-binning by a factor 
of 100 I obtained 263 normal points. The 512 s data comprise 2400 
points of which 291 are near transit. One of the 14 transits was 
rejected due to partial coverage. The remaining 512 s data were 
modelling using JVj nt = 3. A third light of L 3 = 0.07 ± 0.005 
JTinglev et alj201 lh was taken into account. 

The best fits are shown in Fig.[T9]and tabulated in Tables A5 1 
and A52; correlated noise is unimportant. The LD-fit/fix solutions 
of the 32 s data are good, but only the LD-fixed solutions are re- 
liable for the 512 s data. The two ligh t curve solutions agr ee well 
with each other (Table A53) and with lTinglev et alj a201 lh . as do 
the resulting physical properties (Table A54). Improved photome- 
try and spectroscopy would be beneficial. 
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Figure 20. The 5 12 s light curve of CoRoT-15. See Fig.|5]for details. 



6.15 CoRoT-15 

The CoRoT-15 sy stem contains a transiting brown dwarf 
JBouchv et alj|201 II) wit h a mass (65 Mj UD ) ve ry similar to those 
of WASP-30 b (61 Mj up ; lAnderson et alj201 lbh and LHS 6343 Ab 
(63 M Jup ; Sect. l6.27t . CoRoT-15 is the most difficult to study be- 
cause of its faintness (V — 15.5). The CoRoT light curve lasts 
only 31.7 d and covers ten transits at 512 s cadence. Of the 395 dat- 
apoints near transit, four are rejected by a 3<r cut a nd the rest are 
mode lled using JV int = 3 and L 3 = 0.019 ± 0.003 JBouchv et all 

HoI3). 

The JTKEBOP solutions point to a central transit, and therefore 
asymmetric errorbars. The LD-fixed and LD-fit/fix solutions agree 
well (LD-fitted solutions were not tried) so the latter are adopted. 
Unusually for the CoRoT data, correlated noise was found to be 
important, with the residual-permutation errorbars slightly larger 
than the Monte-Carlo ones. The best fit is in Fig.|20]and the model 
details are in Table A55. M y photometric param eters (Table A56) 
are concordant with those of lBouchv et alj J2011I) . I find a similarly 
good match for the physical properties of CoRoT-15 (Table A57), 
which could be improved by further observations of all relevant 
types. 



6.16 HAT-P-4 

We now leave the CoRoT objects behind and turn to a TEP system 
with a light curve from the EPOCH project perfo rmed by the NASA 
Deep Impact satellite ( Christians en et alj|201 lh . HAT-P-4 was dis- 
covered by Kovacs et al. ( 2007) to be a low-density TEP orbiting a 
metal-rich late-F star. lKovacs et alj obtained follow-up photometry 
using the FLWO 1.2 m t elescope and Kepler Cam, which was sub- 
sequently reanalysed by Torres et al. (2008). Further light curves 
were collected bv lWinn et alj ( 1201 lh . who also presented RV obser- 
vations consistent with alignment of the pla netary orbital and stel- 
lar rotational axes. lChristiansen et alj ( 1201 lh analysed the EPOCH 
data, which cover seven consecutive transits followed by three more 

taken five months later. 

The light curve from lKovacs et al . (2007) contains 985 data- 
points spread over two transits, whose uncertainties I have multi- 
plied by 1.67 to get xi ~ 1- The LD-fit/fix solutions are clearly 



better than LD-fixed (larger residuals) and LD-fitted (unphysical 
LDCs). Correlated noise is slightly important. The parameters are 
given in Table A5 8 and the best fit is shown in Fig.1211 

I Winn et alj 120110 observed two transits in the i-band, a full 
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Figure 21. Phased light curves of HAT-P-4 compared to the best fits found 
using JKTEBOP and the quadrat ic LP law. The tra nsit light curves are, 
from top to bottom , z-band from lKovacs et alj 420 07), FLWO i-band from 
IWinn et alj feOllh. FTN t-band fromlwinn et all <20lih . and the EPOCH 
dataset from Christiansen et al. 12011). The residuals are plotted at the base 
of the figure, offset from zero. 



one with FTN and a partial one with KeplerCam. I binned the 662 
points in the first dataset by a factor of 5 to get 133 normal points. 
The second data set was modelled using constraints on the time of 
transit midpoint ( Southworth et alj|2007ah . In both cases LD-fitted 
solutions were not attempted (Tables A59 and A60). 

I rejected one of the EPOCH transits due to insufficient ob- 
servational coverage. Of the 6704 points near the other transits, 68 
were lost to a 4a clip. The data contain substantial correlated noise 
due to pointing wander, so I did not phase-bin them. The residual- 
permutation errorbars are indeed 50% larger than the Monte-Carlo 
ones. Table A61 shows the results, of which the LD-fit/fix were 
adopted. 

The four final light curve solutions are in good agreement 
ixi < 0.04) except for k (Xv — 1-3). and were combined by 
multiplying their probability density functions to get final photo- 
metric parameters (Table A62). These are comparatively imprecise 
because the transit is central, and accord well with literature results. 

The JKTABSDIM results show a significant systematic error 
from different theoretical model sets, but the final results are rea- 
sonable and agree with literature values (Table A63). HAT-P-4 
would benefit from a better light curve. 



6.17 HAT-P-7 



Discovered by |Pal et alj d2008l) . HAT-P-7 was the second TEP 
found in the Kepler field. It is a relatively massive planet (1.8Mj up ) 
orbiting a relatively massive star (1.5 Mg), and was the second 
planet found to h ave a retrograde orbit from RM obse rvations 
JWinn et alj|2009d) after WASP-17 JAnderson et ai1l2010l) . Kepler 
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Figure 22. Phased light curves of HAT-P-7 from the Kepler satellite (upper), 
IWinn et alj (2009a) (middle) and EPOCH (lower) compared to the best fits 
found using JKTEBOP and the quadratic LD law. The residuals are plotted 
at the base of the figure, offset from zero. 



has been observing it at short cadence since the start of the mission, 
and here I have analysed the public data from Q0, Ql and Q2. These 
data have already been used to establish the optical phase curve of 
the pla net (Borucki et al. 2009) and the ellipsoidal effect of the star 
dWelsh et alJuOlO ). The Kepler data have also be subjected to an 
asteroseismic investigation by C hristensen-D alsgaard et al. (2010), 
who used the oscillation spectrum of the data to determine pa to 
high precision. This constraint was not made use of in the current 
analysis, in order to retain homogeneity of approach. The follow-up 
light curve from the discovery paper covers only part of the tran- 
sit and cannot compete with the Kepler data, so was not modelled 
here. 

The Kepler observations are stunning (Fig.l22t and cover 59 
transits (five in Q0, 15 in Ql, 39 in Q2). 37 611 of the 186786 
original datapoints are near a transit, of which 55 were rejected by 
a 4<7 clip and the rest were phase-binned by a factor of 50 to get 753 
normal points. The LD-fitted solutions are good (Table A64) and 
show that both LDCs need to be fitted to account for data with this 
level of precision, although correlated noise is formally important. 

In order to provide a consistency check on the Kepler data 
I have also modelled the ^-ban d observations of one transit (Ta- 
ble A65) bv lWinn et all J2009ah and the EPOCH space-ba sed light 
curve (Table A66) presented by IChristiansen et alj q2010f) . Whilst 
the latter again have substantial red noise, they agree well with the 
Kepler data and the two solutions are therefore combined to ob- 



photometnc p ai 
f Welsh et al.F Ji 



the results of IWelsh et alj ( 120101) is not as good as expected (Ta- 
ble A67). Despite this, HAT-P-7 is one of the best-measured TEP 
systems, alongside HD 209458 (Paper I) andTrES-2 (Sect. l6~28l l. 

The resulting physical properties are contained in Table A68 
and are highly unusual in that the radii of the star and planet 
are measured to precisions approaching 1%. The agreement 
with literature results is g ood o verall but with one caveat. 
IChristensen-Dalsgaard et alj ( 120101) obtained a measurement of 
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Figure 23. Individual transits of HAT-P- 1 1 as seen by Kepler. Successive 
transits are offset by —0.03 flux units for clarity. 



pA = 0.1926 ±0.0023 po by analysing the oscillation spectrum of 
HAT-P-7 A, which is 2.9<r adrift from my transit-derived value of 
PA = 0.2023 ± 0.0024 p© . This discrepancy may either gradually 
sort itself out as additional data appear from Kepler, or alternatively 
might signal the need to include more sophisticated physics in one 
or both of the analysis methods. 

6.18 HAT-P-11 



The third TEP discovered in the Kepler field jBakos et al.1120101) is 
a low-mass planet (0.084 Mj up ) in an eccentric orbit around a low- 
mass star (0.81 Mq). It has been found to have a very oblique or - 
bit through RM observations dWinn et alj2010l:lHirano et al.l20TTh . 
Kepler has observed it in shor t cadence from the start of the mis- 
sion. As with HAT-P-7 above, IChristensen-Dalsgaard et alj a2010h 
analysed the Q0 and Ql data and measured a pa value which has a 
very a high precision but an undetermined accuracy. 

Analysis of the Kepler data is not straightforward, because 
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Figure 24. Phased light curves of HA T-P-11 from the K epler satellite (up- 
per) and in the z-band (lower) from Bak os et alj 120101) . See Fig. 1221 for 
further details. 



there is clear evidence of spot activity on the star. Of the 25 transits 
observed (three in Q0, six in Ql and 16 in Q2), many are marvel- 
lous examples of the phenomenon of a planet transiting a starspot or 
starspot complex and all are affected to some degree (see Fig.|23ll. 
The ratio of the planetary orbital to stellar spin periods is close to 
6.0, so every sixth transit will cross over nearly the same part of the 
stellar surface and thus be similarly affected by spots which evolve 
on a timescale of 29 d or less, a situation that was predicted by 
I Winn et alj ( 120101) . Fig.[23] certainly shows that there are two pre- 
ferred orbital phases for spot activity in the Kepler data. A detailed 
analysis of the spot characteristics will be very interesting. 

In the current work I model the Kepler data using my usual ap- 
proach. The effects of the starspot will therefore be treated as corre- 
lated noise. The oblique orbit of HAT-P- 1 1 b means that it transits, 
at some point, a much greater fraction of the stellar surface than 
an aligned system would, and with a wide distribution of latitudes. 
The assumption that the transited parts of the stellar surface on av- 
erage behave the same as the non-transited areas is therefore more 
justifiable than in cases such as CoRoT-2. Due to LD, the effects 
of starspo ts are much stronger near the centre of the star than its 
limb (e.g. ISanchis-Oieda et alj|201lf) . In contrast, starspot devia- 
tions have the greatest effect on a photometric model when they 
occur at the partial phases of the transit. Treating them as corre- 
lated noise is therefore more likely to lead to an overestimate rather 
than an underestimate of the errorbars. 

The Kepler Q0 to Q2 data comprise 186 802 datapoints of 
which 12 535 are adjacent to a transit. They were not phase-binned, 
as this might affect the estimation of the starspot-induced cor- 
related residuals of the fit. The Kepler data don't agree exactly 
with previous orb ital ephemerides of the system so the To from 
iBakos et alj q2010 | > was included as a constraint following the ap- 
proach of ISouthworth et alj l2007af) . I also accounted for orbital 
eccentr icity using e cos cj = 0.261 ± 0.082 and e sin uj = 0.085 ± 
0.043 l lWinnetal.ll2010h . The LD-fitted solutions are poor and 
have unphysical LDCs so the LD-fit/fix solutions were preferred 
(Table A69). The residual-permutation errors were larger than the 
Monte Carlo ones for the inclination but not for the other photo- 
metric parameters, supporting the approach taken to fit the data. 
The best fit is in Fig.1241 
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Additional light curves are available from lBakos et alj (2010). 
comprising eight transits observed in the 2 -band using KeplerCam, 
one in the r-band with the same instrument, and three in the /-band 
with the Konkoly Schmidt. The first of these datasets is worth solv- 
ing as a check of the Kepler data: the 4110 datapoints were phase- 
binned by a factor of ten to obtain 410 normal points. Correlated 
noise was not found to be important, and the LD-fit/fix solutions 
were adopted (Table A70). The agreement with the Kepler data is 
good (Table A7 1). Literature results are in moderate agreement, but 
the current results are to be preferred as they are the first to be based 
on the Kepler observations. 

The JKTABSDIM results suggest that the measured T e s of 
the host star is quite high (4780 ± 50 K), causing the system to 
occupy the zero-age edge of the model grids. A hefty decrease 
of 300 K in T c g would be needed to assuage this problem, and 
would shift Ma from 0.82 to 0.75 Mq with other parameters less 
strongly affected. I therefore acquired some alternativ e r e g esti- 
mates from a number of sources. The Tycho 73— y value dH0g et al.l 
1 19971) supplemented by the calibration of ISousa et al.l ( 120081) re- 
turns T e g — 4852 K. B — V is not generally regarded as a good 
T e g indicator for late- type dwarfs: the same value for this colour 
index equates to a K5 star with T e g = 4410 K using the tables of 
IZombeckl Jl990h . Using the 2MASS JHK S dSkrutskie et al.l2006h 
magni tudes is a better be t : the V — K s colour index and the calibra- 
tion of lCasagrande et al.ld2008ri yields T cff = 4765 K. Dr. B. Smal- 
ley has kindly calculated the luminosity of HAT-P- 1 1 A from the 
available broad- band optical and infra red photometry and the Hip- 
varcos parallax jPerrvman et al.ll997i) to be logfL/ Lq) = —0.61. 
This agrees with the JKTABSDIM solutions for the measured T B g of 
the star (-0.63) but not with those for the l ower T cff (-0.77). The 
T c ff measurement from lBakos et alj ( 120101) is supported by the in- 
vestigations above. 

Table A72 shows my calculated physi cal properties fo r HAT- 
P-llb compared to literature solutions. Bakos et al.l d20101) find 
smaller radii for the two components which may be down to dif- 
ferences in analysis. In the HAT methodology the stellar param- 
eters are forced to agree with a theoretical stellar model (usu- 
ally interpolated from the Y 2 isochrones); in my solution process 
I find the point of closest agreement but do not require this to ex- 
actly reproduce a point in a grid of theoretical predictions. I find 
a very young but poorly constrained age, which is in accordance 
with the starspot activity seen in both the Kepler and the HAT 
data. IChristensen-Dalsgaard et al . (2010) derived pa = 1.7846 ± 
0.0006 po (errorbar does not account for systematic errors) from 
the oscillation spectrum of HAT-P- 11 A in the Kepler light curve. 
This indicates a much less dense star than I find (2.415±0.097 Pq), 
a similar situation to HAT-P-7 but with a much stronger discrep- 
ancy. Further investigation is needed to understand the disagree- 
ment; extra RV measurements would be useful to refine the K\, 
e cos uj and e sin ui values. 

6.19 HD 17156 

Th e planet orbiting HP 17156 was discovered by the RV method 
by iFischer et alj d2007D . who also searched for transits but did 
not detect them or rule out their existence. A transit w as soon 
observed using three telescopes by Barbie ri et alj J20071). as par t 



2003) 



2008) 



of the transitsearch.org network (Seagroves et al. 
Further transi t light c urves were obtain e d by [Gillon et alj i 
iBarbieri et ail J2009I) an d llrwin etaD d2008h . Two good light 
curves were presented by I Winn et alj d2009ch . covering the tran- 
sit on Christmas Day in 2007. The RM effect has been observed 
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Figure 25. The light curves of HD 17156 compared to the JKTE BOPbest fit 
Top is H ST/FGS, middle and bottom are (6 + y)/2 and z from lwinn et alj 
d2009cb . See Fig.Fj2]for further details. 



on multiple occasions (Naritaetal. 2008; Cochran et al. 2008; 
Barbieri et al. 2009; Narita et al. 2009) and overall indicates axial 
alignment. 

An extensive a naly sis of HD 17156 w as performed by 
iNutzman et"all d201 lh and lGilliland et alj d201lh ba sed on data ob- 
tained using the HST/FGS. The asteroseismic study dGilliland et alj 
2011) yielded a mean density and age for the star. The theoretical 
uncertainty of the former quantity is probably small, but is signifi- 
cant for the latter q uantity. The HST/FGS data cover three transits, 
which were used bv lNutzman et alj d201ll) to measure the physical 
properties of the system. Analyses including the asteroseismic p\ 
as a constraint are in good agreement with standard analyses, but 
with smaller errorbars (by a factor of three for R\ and R^). I did 
not include the asteroseismic constraint in my analysis, in order to 
retain homogeneity with the results for other s ystems. 

I have solved the two light curves from Iwinn et alj ((2009c), 
which are in the (b + j/)/2 and z passba nds, each cover i ng on e 
transit. I also modelled the HST data from lNutzman et alj ( 1201 lh . 
which comprises partial coverage of each of three transits, using 
a light curve which has been corrected for instrumental effects. 
The FGS data were taken with the F583W filter which covers 440- 
710nm, so I adopted for LDCs appropriate for a combined g+r fil- 
ter. The orbital shape of HD 17156 was accounted for by adopting 
e = 0.6768 ± 0.0304 and u = 121.71 ± 0.43 deg dNutzman et ail 
l201ll) . The best fits are exhibited in Fig.l25l 

The results for the HST light curve are given in Table A73 and 
LD-fit/fix is the best alternative. Correlated noise is marginally im- 
portant. The (6 + y)/2 data prefer a higher i and thus smaller r\ 
and »*b (Table A74), whereas the z observations yield the opposite 
situation (Table A75). The two datasets agree with the HST results 
overall, so the HST ones are adopted as the final parameters (Ta- 
ble A76). Their agreement with published values is reasonable. 

The physical properties of HD 17156 are shown in Table A77 
and reveal some model-dependent error which stems primarily 
from the DSEP models. The agreement with literature values is 
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Figure 26. The Spitzer light curve of HD 80606 (Hebrard et al. 20KJ). See 
Fig.|22]for further details. 

again good, and I provide the first measurement of T B ' q and O. 
HD 17156 is a well-characterised system which is not in need of 
further observations. 



6.20 HD 80606 

HD 80606 b was found to be a p lanetary-mass object with a period 
of 1 11.8 d by iNaef et alj ( 1200 II) from extensive RV observations. 
At the time o f discovery it w as the most eccentric exoplanet known 
(e = 0.927). iLaughlin et al. (2009) observed an occultation of the 
planet whilst using Spitzer to probe the heating of th e planet around 
periastron. Based on the parameters then known, ILaughlin et all 
calculated a probability of 15% that a transit also occurs, and en- 
couraged follow-up observations to detect it. The system was mon- 
itored over the next predicted ti me of transit and duly found to 
be a T EP by Moutou et al. ( 2009). lGarcia-Melendo & McCulloughl 
J2009h and lFossev etal] J2009h . 

RM observa t ions of HD 80606 have be e n obtained and s tudied 
by iMoutou et all d2009t> . iPont et all J2009I) . IHebrard et ail J2010h 
and I Winn et all ([2009b), with a good agreement that the orbit is 
oblique to high confidence. Coupled with the high ecce ntricity, this 
has b e aring on its formation and evolution mechan isms dNaef et all 
1200 ll ; lMardling&Linll2002l; iMatsuo et all |2007|) . The relatively 
high mass of the planet (4.11 Mj up ) agrees with the established 
pattern that TEPs in ecce ntric orbits are generally the massive ones 
JSouthworth et alll2009cl) . 

Spectral sy nthesis ana ly ses o f the parent st ar have been 
performed by iNaefetalJ J200lh . ISantos et all J2004) and 
iGonzalez et"ail J2010I) : I adopt the T eB from ISantos et alj J2004) 
as it is the middle of th e three values and its e rrorbar encompasses 
the other two; I use the lGonzalez et alj J2010I) [^1 value because 
all three agree and it has an equitable errorbar. I take the Ka value 
from the most recent detailed RV analysis, Winn et al. (2009b). 

Good light curves of HD 80606 are difficult to obtain 
due to the long transit duration (nearly 12 hours) and or- 
bital period. The transit discovery l i ght c u rves (M outou et al. 
l2009l : iGarcia-Melendo & McCulloughl boosl : iFossev et alj 120091) 
all mis sed ingress as it occurred during daylight. Wi nn et all 
(2009b) combined data from eight observing sites spread through 
mainland USA and Hawaii to obtain full coverage of the 2009 June 



tran sit, but the resulting l ight curve was still a lon g way from defini- 
tive. lShporer et al.l ( l201CJ) and Hidas et alJJ2010l) performed similar 
observing campaigns covering several transits, but the data are un- 



fortunately strongly affected by systematics. IHebrard et all 120101) 
used Spitzer to obtain a complete and high-quality light curve of 
the 2010 January transit in the IRAC 4.5 /im passband. This is by 
far the best light curve of any transit of HD 80606, and is the only 
one analysed here. 

The Spitzer data were binned by a factor of 100 to lower the 
number of datapoint^j from 31767 to 318. The resulting time 
resolution is 215 s, which is reasonable for such a long-duration 
transit. A disturbance near mid-eclipse is noticable and can be at- 
tributed to starspots (Fig.l26t; this was treated as correlated noise 
in the JKTEBOP analysis (see Sect. l6.18t . The flux from the planet 
in this near-infrared passband must be accounted for, which was 
done by incl uding a light ratio of 0.00 1 ± 0.001 with the method 
described bv lSouthworth et alj d2007al) . This light r atio gives a sec- 
ondary eclipse of similar depth to that observed by Laughlin et al. 
(2009), with a conservative errorbar to allow for inefficient en- 
ergy redistribution within the planetary atmosphere. The orbital 
shape was accounted for usin g ecosu} = 0.4774 ± 0.0018 and 
esmui = -0.8016 ± 0.0017 JHebrard et alfcoiol) . The LD-fitted 
solutions yield unphysical LDCs so the LD-fit/fix solution was 
adopted (Table A78). Correlated noise is important, due to the treat- 
ment of the starspot. The final photometric parameters (Table A79) 
are in reasonable agreement with published values. 

The JKTABSDIM results show that the systematic error (model 
disagreement) is relatively high, and dominates the error budget for 
Ma, Mb and a (Table A80). The solution with the dEB constraint 
is also somewhat different to the solutions using theoretical models. 
The source of these problems is not obvious, but the discrepancies 
are small enough that HD 80606 is one of the best-characterised 
TEP systems. 



6.21 Kepler-4 

Kepler-4 was the first TEP system to be announced as a discovery 
of the Kepler satellite (Borucki et al. 2010a) and is a very low-mass 
object (0.075 Mj up ) orbiting a slightly evolve d star (l og g = 4.17). 
It has also been studied by Kipping & Bakosl 1201 lbl) . Eccentricity 
is significant at only 2a and additional RV measurements are nec- 
essary to investigate this further. The available Kepler data cover 
13 shallow transits at long cadence (29.4 min) in Q0 and Ql and 26 
transits at short cadence (58.8 s) in Q2. 

The addition of Q2 data over Q0 and Ql allows an improve- 
ment of the orbital ephemeris. The short-cadence data were binned 
to the time resolution of the long-cadence data. The long-cadence 
and binned short-cadence data were then modelled using JKTEBOP 
and iVint = 10 to find: 

To = HJD 2 454 956.61132(92) + 3.213658(38) x E 

where the bracketed quantities give the uncertainties in the preced- 
ing digit, and the uncertainties come from 100 MC simulations. 

The long-cadence transit data were cut from the full light 
curve and normalised as described in SectJJJ] giving 323 datapoints 
out of the original 2101. The 123 536 short-cadence observations 
were treated similarly, then phase-binned by a factor of 25 to yield 
708 normal points. I assumed a circular orbit and a third light of 
I/3 = 0.02 ± 0.02. The long-cadence data do not constrain the fit 
well; parameter perturbations could not be applied in the Monte 



14 The Spitzer data I used had already been binned by a factor of five before 
being lodged with the CDS. 
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Figure 27. Phased long-cadence (top) and short-cadence (second from top) 
light curves of Kepler-4 compared to the best fit found using JKTEBOP and 
the quadratic LD law. For the long-cadence data the fit is shown by a blue 
line, and an evaluation of the same model but without numerical integration 
is shown by a purple line. The residuals are plotted at the base of the figure, 
offset from unity. The purple line through the residuals shows the difference 
between the fit to the data and an evaluation of the same model but without 
numerical integration. 



Carlo analysis so these errorbars may be optimistic. The long- 
cadence data were solved using N[ nt — 10 (Table A81) and the 
short-cadence data were solved without using numerical integra- 
tion (Table A82). The best fits are given in Fig.[27] In both cases the 
LD-fit/fix solutions were best. The residual-permutation errorbars 
were much larger than the Monte-Carlo ones for the long-cadence 
light curve, which is likely due to the inability to apply parameter 
perturbations in this case. 

The solutions for the two datasets agree well and were com- 
bined by multiplying their probability density functions to give fi- 
nal photom etric parameters (Table A8 3). The photometic param- 
eters from Kipping & Bakosl ( 1201 lbf) disagree with those from 
iBorucki et al.l C2010ah : the latter find a low er inclination and thus 
larger ta and rv My results support those of Borucki et al. (2010a) 
and I a m unable to completely reproduce the Kippi ng & Bakosl 
( 1201 lbl) results by modifying my treatment of eccentricity or nu- 
merical integration. 

The physical properties of Kepler-4 (Table A84) show a rea- 
sonable agreement with literature results albeit with substantially 
larger errorbars in some cases. This is despite the availability of 
much more extensive (and higher-cadence) observations available 
for my analysis. The Kepler satellite continues to observe the sys- 
tem at short cadence so a much improved light curve will gradually 
accumulate. Additional spectroscopy to improve the T e g and Ka 
measurements should be a high priority. 



6.22 Kepler-5 

The discovery of Kepler-5 was announced by iKoch et al.l 020100 : 
the planet is relatively massive (2.04 Mj up ) and the star is 
quite evolved (logff = 4.17). It has also been studied by 
iKipping & Bakosl J2011br) . The available Kepler data now comprise 
12 transits studied at long cadence during Q0 and Ql, of which 330 
of the 2101 datapoints occur near to a transit, and 23 transits at 
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Figure 28. Phased light curves of Kepler-5. See Fig.|27]for further details. 



short cadence after rejection of one transit due to systematic noise. 
The 123 536 short-cadence datapoints were reduced into 379 nor- 
mal points. A revised orbital ephemeris was measured from the Q0, 
Ql and Q2 data: 

T = HJD 2 454 955.90059(36) + 3.548469(15) x E 

using the same approach as for Kepler-4. 

A th ird light of L3 = 0.02 ± 0.002 and a circular orbit were 
assumed (Koch et al. 2010) and the long-cadence data were mod- 
elled using 7V"i n t = 10. The JKTEBOP solutions (Tables A85 and 
A86) show that correlated noise is unimportant and that the LD- 
fitted results are viable. Agreement between the two datasets is not 
good and can be attributed to information loss at the long cadence. 
I therefore adopt the short-cadence solutio ns as final (Table A87). 
I find a notably s maller ta and r^ than IKoch et alj d2010l) and 
IKipping & Bakosl J2011bl) . which reflects the difference between 
the long-cadence data (which were available for those studies) and 
the better short-c adence obse r vation s (which were not). The uncer- 
tainties given by IKoch et al.l ( 1201 01 are too low. The best fits are 
shown in Fig.1281 

The JKTABSDIM results are contained in Table A88 and ex- 
pectedly disagree with literature studies based on only the long- 
cadence data: I find smaller masses and ra dii for both star and 
planet. The discovery paper dKoch et alJuOl O) obtained two solu- 
tions corresponding to different evolutionary stages for the host star 
(with masses 1.21 and 1.38 Mq), and endorsed the more evolved 
alternative. I do not find this problem, due to my refined photomet- 
ric parameters. Additional photometric observations are needed, 
are currently being obtained by the satellite, and will make Kepler- 
5 one of the best-charaterised TEP systems. 



6.23 Kepler-6 

Kepler-6 was found to be a TEP by iDunham et al.l J20101) and 
is interesting because of the high metallicity of the host star 
/[if] = 0-34±0.05). The same long-cadence data were studied by 
IKipping & Bakosl J2011bf) , who found a solution with a higher or- 
bital inclination and 2.6cr signficant signal in the TTV peridogram 
which may indicate stellar activity. The Kepler Q0 and Ql data 
cover 13 transits at long cadence, whereas the Q2 data include 41 
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Figure 29. Phased light curves of Kepler-6. See Fig. l27l for further details. 



transits at short cadence. A new orbital ephemeris was determined 
from these data: 

To = HJD 2 454 954.485805(64) + 3.2347020(33) x E 

The long cadence data comprise 2102 points of which 246 
are near a transit. These were modelled using N-mi, = 10 to yield 
the results in Table A89. The solutions are very sensitive to the 
treatment of LD: LD-fixed gives a poor internal agreement and 
LD-fitted returns unphysical LDCs, so the LD-fit/fix results were 
adopted. The 123 536 short-cadence datapoints were reduced into 
34 normal points and the LD-fit/fix so lutions were best (Tab le A90). 
In both cases L s = 0.033 ± 0.004 JDunhamet al]|2010h was in- 
corporated and correlated noise was found to be inconsequential. 
The two datasets agree well so the parameters were combined (Ta- 
ble A91). Published studies agree well with my results, although 
their errorbars are questionable. The best fits are shown in Fig.|29l 

The physical properties of Kepler-6 are given in Table A92. 
The DSEP model solutions disagree with the others so are not in- 
cluded in the fin al results. I find a so mewhat s maller planet and 
star com pared to lDunham et al.l (2010), whereas lKipping & Bakosl 
( 120110) agree with my results within the errors. Kepler continues to 
observe Kepler-6 at short cadence and additional RV measurements 
would also be useful. 



6.24 Kepler-7 

Discovered by lLatham et al.l 12010D . this a very low-density TEP 
(0. lOpjup) around a slightly ev olved star (lofiff = 3.96 ). The 
same data were also studied by Kip ping & Bakosl (201 lbl) , who 
detected an occultation with a significance level of 3.5<r. The Q2 
data are long-cadence, so no short-cadence data are available for 
analysis. The Q0, Ql and Q2 observations cover 26 transits, and 
947 of the 6177 datapoints were retained for the JKTEBOP analysis. 
1/3 = 0.025 ± 0.005, a circular orbit and Mnt = 10 were adopted. 
The updated orbital ephemeris is: 

To = HJD 2 454 967.27598(11) + 4.8854948(82) x E 

The LD-fitted results are reliable and have reduced residuals com- 
pared to the LD-fit/fix solutions (Table A93). Correlated noise is 
insignificant. I find a lower inclination and therefore a higher r\ 
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Figure 30. Phased long-cadence light curve of Kepler-7. See Fig. l27l for 
further details. 
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Figure 31. Phased light curves of Kepler-8. See Fig.|27]for further details. 



and rb than previous studies (Table A94), and investigations reveal 
that this is primarily due to inclusion of the Q2 data. The best fits 
are shown in Fig.|30| 

The physical properties of Kepler-7 (Table A95) reveal that 
the star is one of the most evolved known to host a TEP, and 
the planet itself is the second-m ost rarefied known after WASP- 
17 (Anderson et aI.ll201Ctl201 laf) . As expected from the photomet- 
ric parameters, my physical properties do not agree well with those 
previously published. Confirmation of this will be possible soon, as 
Kepler-7 has been observed in short cadence by Kepler from Quar- 
ter 3 onwards. New RV measurements would also be worthwhile. 



6.25 Kepler-8 

Like Kepler-7, Kepler-8 is a low-densit y TEP (0.21 p j UD ) o rbiting 
a slightly evolved star (logg = 4.18). I Jenkins et al.l J2010al) and 
Kipping & Bakos (201 lb) have studied the Kepler Q0 and Ql long- 
cadence data, whereas the Q2 short-cadence observations are now 
available. A refined orbital ephemeris was established as: 

To = HJD 2 454 954.11844(18) + 3.5225047(76) x E 

The long-cadence observations cover 13 transits, and I retain 
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Figure 32. Phased light curve of KOI-428. See Fig.|27]for further details. 
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Figure 33. Phased light curve of LHS 6343. See Fig.[2Z]for further details. 



265 of the original 2098 datapoints (solved using iVmt = 10). 
The Q2 data consist of 123 536 points covering 24 transits, which 
I phase- binned down to 299 normal points. The third light value 
given bv ljenkins et al.l fcOlOah does not have an errorbar, so I adopt 
1/3 = 0.0075 ± 0.0075 to be conservative. The best fits can be seen 
in FigJJT] Models of the long-cadence data are rather sensitive to 
the treatment of LD and are also in comparatively poor agreement 
with the short-cadence models (Tables A96 and A97). I therefore 
adopt the LD-fit/fix flavour of the latter (Table A98), which is also 
in accord with and more precise than literature studies. The error- 
bars found by Jenkins et al. (2010a) are smaller than expected. 

The absolute dimensions of Kepler-8 (Table A99) agree well 
with published values and establish it as a well-understood system. 
Kepler continues to observe Kepler-8 photometrically and further 
RV observations are merited. 



6.26 KOI-428 

KOI-428 was one of the 306 Kepler Obje cts of Interest that was pre- 
sented to the astronomical community bv lBorucki et alJa201 lh ; the 
400 brightest ones were retained for in-house anal ysis by the Ke- 
pler T eam. Follow-up spectroscopic observations bv lSanterne et al.l 
( 120 111) have subsequently shown this to be a system containing a 
relatively massive planet (2.1 Mj up ) orbiting a comparatively hot 
star (65 10 K). 

The Kepler Ql and Q2 data contain 12 transits observed at 
long cadence, compa red to only four in the Ql data available to 
ISanterne et alj 1201 ll) . 741 of the 5708 datapoints were retained; 
one of the 12 transits was rejected because it is deformed by instru- 
mental artefacts arising from a pointing jump. These were solved 
using iVint = 10, a circular orbit and no Ls (Table A100). Corre- 
lated noise is negligible. I find a lower inclina tion and thus larger 
r\ and r^ compared to ISanterne et alj ( 1201 lh . but the values are 
within their errorbars (Table A101). Fig.[32]shows the best fit to the 
light curve. The revised orbital ephemeris is: 

T = HJD 2 455 005.51858(50) + 6.873130(75) x E 

JKTABS DIM returns a notica bly larger but less massive star 
compared to lSanterne et alj (1201 lh . as expected given the slightly 
different photometric parameters (Table A102). The model fits pre- 
fer a T e ff lower by 80 K (0.8<r), so this may indicate that analysis 
of the Kepler data in isolation results in a star which is too large. 
KOI-428 remains on the Kepler target list so additional data will be 



available soon. The errorbars quoted bv lSanterne et alj ( 1201 lh are 
too small. Further spectroscopic study of this object is warranted. 

6.27 LHS 6343 



LHS 6343 is a nearby M dwarf JLuvterj 19791) which was found 
to host a small transiting object by Ijohnson et al.l 1201 lh based 
on Ke pler data. Follow-u p imaging and spectroscopic observa- 
tions J Johnson et alj|201 lh revealed that the system contains two 
M dwarfs (A and B) separated by 0.55" and that the brighter com- 
ponent A hosts a likely brown-dwarf companion with an orbital 
period of 12.71 d. The Kepler Q0, Ql and Q2 data contain a total 
of 1 1 transits observed at long cadence. 

The Kepler observations are of the combined flux of the three 
components, and a light ratio of B versus A has not been directly 
observ ed. The adaptive-optics imaging obtained by Johnson et al. 
1201 lh give magnitude differences in the JHK S passbands of 
0.49 ± 0.05, 0.49 ± 0.05 and 0.45 ± 0.05, respectively. The char- 
acteristics of these passbands are such that an extrapolation to the 
wide Kepler passband is reasonable for stars as similar as A and 
B. This was performed usin g the process outlined in PapeiTII (see 
also lSouthworth et alj|20ld) . resulting in L 3 = 0.29 ± 0.09. Mea- 
surements of the flux ratio at optical wavelengths would be useful 
in re fining these numbers. 

Ijohnson et ail d201ll) obtain T eB = 3130 ± 20 K for com- 
ponent A from the calculated V and observed K a magnitudes of 
the sta r applied to the photometric calibrations of lCasagrande et al.l 
(2008). The quoted errorbar is clearly a precision rather than a true 
uncertainty. I therefore extrapolated the infrared light ratios into the 
V band to obtain the V — K a colou r index which, with the calibra- 
tions from lCasagrandeetal.1 J2008h . gives T cB = 3300 ± 200 K 
for component A where the errorbar is conservative. My T e s value 
was used in the analyses below. 

The 6177 long-cadence datapoints from the Kepler Q0, Ql 
and Q2 observations were slimmed down to 268 points nearby a 
transit. These were modelled with JKTEBOP using L3 = 0.29 ± 
0.09 (see above), e = 0.056 ± 0.032, uj = 337° ± 56° 
( Ijohnson et alj201 lh and 7Vi nt = 10. 1 included P or b and Tp as fit- 
ted pa rameters constrained by the To value found bv ljohnson et al.l 
1201 lh from a ground-based observation of one transit. The ensuing 
orbital ephemeris is: 

To =HJD 2 454 995.358014(44) + 12.7138107(73) x E 

The LD-fixed solutions are poor (TableA103). The LD-fit/fix al- 
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Figure 34. Phased light curves of TrES-2 compared to the best fits found 
using JKTEBOP and the quadratic LD law. From top to bottom the data are 
from Kepl er, the 790.2 nm th en 794.4 nm sets fr omlColon et all 120101), R - 
band from Rabus elail <2009l) . and EPOCH from lChristiansen et aljT201lh . 
See Fig. l26l for further details. See Paper III (sect. 4.7 and fig. 11) of Paper III 
for previous results. 



ternatives are better, but the quality of fit is hampered by scat- 
tered data mischievously placed right at the transit midpoint. Ke- 
pler continues to observe LHS 6343 so an improved light curve 
is in the pipel i ne. Co rrelated noise is not important. Compared to 
Ijohnson et al.l (1201 lh I find a smaller k and larger n, as well as 
larger errorbars despite having many more photometric observa- 
tions (Table A104). The best fit is in Fig.|33l 

Derivation of the physical properties of the star and its sub- 
stellar transiting companion is difficult due to the low mass of 
the former object. Theoretical stellar models are unreliable in this 
regime (see fig. 4 in Paper III) so an additional systematic error 
should be ad d ed to those quoted in Table A 105. Compared to 
Ijohnson et al.l J2011I) I find the star to be more massive, which 
propagates into a correspondingly larger mass for the companion of 
Mb = 70 ± 6 Mj U p close to the stellar/substellar boundary. Further 
spectroscopic and spatially resolved optical observations would be 
useful in pinning down the T e fl of the host star and the mass of its 
companion. 

6.28 TrES-2 

TrES-2 was discovered by the Trans-Atlantic Exoplanet Survey 
dO'Donovan et al.1120061) and subsequently treated in Paper I and 
Paper II. TrES-2 was the first TEP discovered in the Kepler field, 
and a l ight curve of stunni ng quality from this satellite is now avail- 
able ( Gilliland et al. 2010). An analysis of these data has been given 
by I Kipping & Bakosl 12011a ). In addition, excellent ground-based 



light curves have been obtained by IColon et al.l ((2010) and it was 
one of the targets of the EPOCH project u sing the NASA Deep 
Impact spacecraft ^Christiansen et al.ll201 ll) . We also add to this 
datase t the Johnson i?-band light curve obtained by iRabus et al.l 
( 2009) in the course of a TTV study. The sys tem offers one compli - 
cation: a fainter star at a separation of 1.09" JDaemgen e t al. 2009). 
TrES-2 was revisited in Paper III to account for this situation. 

The Kepler data cover four transits in Ql, 14 in Ql and 34 
in Q2 (ignoring one near an instrumental artefact), all at short 
cadence. The 186 802 original datapoints were reduced to 18 310 
by rejecting observations far from transit, then to 18 297 by a 4a 
clip, and then to 311 normal points by phase-binning by a fac- 
tor of 30. L 3 = 0.0258 ± 0.0008 was used and the LD-fit/fix 
solutions adopted (Table A106). Correlated noise was not impor- 
tant. Although Kepler continues to observe TrES-2, it is already the 
photometrically best-measured TEP known. The light curve fits are 
shown in Fig.[34] which is best viewed in conjunction with Fig. 1 1 
in Paper III. 

The two light curves from lColon et al.l J20101) cover the same 
transit in two very narrow passbands, at 790.2 and 794.4 nm, ob- 
tained using the 10.4 m GranTeCan and OSIRIS imager equipped 
with a tunable filter. L 3 = 0.0355 ± 0.0005 was used for both 
passbands. For the 790.2 nm data the LD-fixed solutions had to be 
adopted, and correlated noise was found to be moderately impor- 
tant (Table A 107). The LD-fit/fix solutions could be adopted for the 
794.4 nm data, for which correlated noise was found to be unimpor- 
tant (Table A 1 08 . 

The i?-band data from lRabus et al . (2009) comprise data from 
five transits which are supplied already phased and binned. I scaled 
the errorbars by a factor of 0.485 to obtain ^«1 and modelled 
them using L 3 = 0.0287 ± 0.0007. The LD-fit/fix solutions are 
reasonable (Table A 109) and correlated noise is incidental (as ex- 
pected given the phase-binning process). 

The EPOCH data (Christian sen et aLlPJoTlh cover eight tran- 
sits. 3534 of the original 27 724 datapoints were retained, of which 
17 were subsequently rejected by a 3.5<r clip. They were not phase- 
binned as correlated noise was correctly expected to be important. 
1/3 = 0.026 ± 0.002 was assumed. I had to adopt the LD-fixed so- 
lutions as the LD-fit/fix and LD-fitted alternatives gave anomalous 
results and negative LDCs (Table Al 10). 

The overall results for each light curve are given in Table A 1 1 1 
and show an excellent agreement overall, k is as expected the least 
concordant parameter, but even here the agreement is at the level 
of Xv = 0.9. The final photometric parameters are the weighted 
means of the ones for each light curve. Their agreement with pub- 
li shed values is in general e xcellent except, perplexingly, for those 
o flKipping & Bakosl l2011ah which are the only other ones to be de- 
rived from the Kepler data of TrES-2. Similar concerns have been 
noted for Kepler-4 to Kepler-8, so there may be a small systematic 
difference in the results from Kipping & Bakos compared to other 
researchers. 

TrES- 2 is of par t icular interest because iMislis & Schmittl 
(2009) and lMislis et alj yOlOJ) have found evidence for a decrease 
in the system's orbital inclination. The Kepler data rule out an ef- 
fect of the expected size, and do not provide ev idence of changes 
in any of the photometric parameters (see also Kipping & Bakos 
l2011al) . A natural explanation of the previous detection of a change 
in inclination would be the presence of subtle systematic errors in 
transit light curves. 

The JKTABSDIM results are given in Table All 2 and show that 
TrES-2 is now very well characterised. The photometric parame- 
ters contribute only a small part of the error budget for the mea- 
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Figure 35. Phased light curves of TrES-3 compared to the best fits found 
using JKTEBOP and the quadratic L P law. From top t o bottom are the 
790.2 nm then 794.4 nm datasets fro m lColon et a?] feOlfJ) . and the EPOCH 
data from lChristiansen et all feplll) . See Paper III (sect. 4.8 and fig. 12) for 
previous results. See Fig.[26]for further details. 



surements of its physical properties. The best way to improve the 
results further would be to obtain a more precise M|m value. 



6.29 TrES-3 



TrES-3 was discovered bv lO'Donovan et al.l ( 120070 and previously 
studied in Paper III. Since then interleaved light curves of one tran- 
sit have been published in two narrow passbands, with central 
wavelengths 790.2 nm and 794.4 nm, bv JColon etafl (J2 010). pho- 
tometry of four transits has been presented bv lLee et alJd201 11) . and 
a light curve from the Deep Impa ct mission has been obta ined in 
the course of the EPOCH project dChristiansen et alj|201 11) . In the 
current work I improve upon the results f rom Paper III (see fig. 12 
in that work) by s tudyin g the data from IColon et alj ( 120101) and 
IChristiansen et alj j20lTh . 

The 790.2 nm and 794.4 nm light curves cover the same transit 
by alternating between the two passbands using a tunable filter. In 
both cases (Tables Al 13 and Al 14) correlated noise is unimportant 
and LD-fit/fix provides the best solution. 

The EPOCH data cover six transits, of which one has only 
partial coverage and one has problems with systematic noise. The 
1171 datapoints in the region of the remaining four transits were 
selected and four were rejected by a 3.5a clip. Table Al 15 shows 
that the LD-fit/fix solutions return questionable results, so the LD- 
fixed solutions had to be adopted. Unusually for the EPOCH data, 
the residual-permutation errorbars were not larger than the Monte- 
Carlo ones. 

The final results for each light curve are given in Table Al 16, 
alongside the results for the seven light curves investigated in Pa- 
per III. Fig.[35]shows the best fits. The final photometric parameters 
were obtained by multiplying the probability density functions for 
the ten light curves from Paper III and the current work. The agree- 
ment between light curves and with literature values is excellent. 
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Figure 36. Phased light curves of WASP-3 compared to the best fits found 
using JKTEBOP and the qua dratic LP law. From t op to bottom the datasets 
are q then i then z from | Tripathi et alj J2010h . and the EPOCH data 
from lChristiansen et alj 1201 II) . See Fig. l26l for further details. See Paper III 
(fig. 14 and sect. 4.10) for previous results. 



The physical properties of TrES-3 can be found in Table Al 17 and 
again show good correspondance with published numbers. TrES-3 
is well-characterised, but would benefit from a few more RV mea- 
surements. 



6.30 WASP-3 



WASP-3 was treated in Paper III but is revisited here, as i t has 
since been obser ved from th e ground bv| Tripathi et al. (2010) and 
from space by IChristiansen et al. (2011 ). Six nice transit light 



curves were obtained by iTripathi et alj J2010I) . comprising three 



in i and two in g from the FLW O 1.2 m and one in z fro m the 
University of Hawaii 2.2 m. The IChristiansen et alj ( 1201 lh data 
come from EPOCH and cover eight transits with good precision 
but som e systematics due to poin ting wander and imperfect flat- 
fielding. iMacieiewski et alj ( 120 10 ") has detected possible TTVs in 
the WASP-3 system, which are yet to be independe ntly confirmed. 

I adopt the same T e s and [Epl as in Paper I dPollacco et al.1 

120081) . Two RV studies exist f or WASP-3, yielding velocity am- 
plitudes of 290.5 ± 9.5 ms" 1 JTripathi et al]|20ld) and 278.2 ± 
13.6 ms" 1 dMilleret alJl2O10h . In Paper III I took the former of 
these two values, but in the current work I adopted instead the 
weighted mean: Kx - 286.5 ± 7.8 ms" 1 . 

The gr-band data contain two transits so the P or b was included 
in the fit to insure against possible bias from TTV effects. The LD- 
fit/fix solutions are best and the residual-permutation errorbars are 
about 30% larger than the Monte-Carlo ones (Table Al 18). The i- 
band data cover three transits so P or b was again fitted for. The LD- 
fit/fix solutions are good and correlated noise is not important (Ta- 
ble Al 19). The z-band data encompass one transit; LD-fit/fix val- 
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ues were adopted and the residual-permutation errorbars are 25% 
larger than Monte-Carlo (Table A120. The EPOCH data cover eight 
transits, of which one was rejected due to poor observational cov- 
erage, with 18 622 datapoints. 3397 are near transit, of which 34 
fell foul of a 3.5ct clip. The LD-fit/fix results were retained and the 
residual-permutation errorbars are 40% larger than the Monte Carlo 
ones (Table A121). 

Fig.1361 shows best fits of the four light curves. Table A122 
summarises the photometric results from the current work and from 
Paper III. All seven light curves are used to calculate the final pho- 
tometric parameters, which are mostly in good agreement. The k 
values are more scattered than they should be (\u — 5-4) and this 
has a knock-on effect on rb which has been accounted for in the er- 
rorbars. Other works are in good agreement with my results, albeit 
with unreasonably small errorbars in some cases. An exception is 
iMiller et al.l (120100 who find a discrepant solution with high incli- 
nation and thus lower rx- This situation propagates into the physi- 
cal properties (Table A123), which are now substantially improved 
over those in Paper III. Further spectroscopic observations, both for 
RV and atmospheric parameter measurements, are warranted. 

6.31 Other TEPs 

I have returned to the XO-4 system, which has received a new and 
substantially improved Ka measurement from lNarita et al.l ( 120 id) 
since Paper III. Table A 124 shows the revision in the system param- 
eters this brings. 

Finally, I have checked that the modifications to the JK- 
TABSDIM code outlined in Sect. |3.1| (primarily the much better 
sampling in age) by re running solutions for the WASP-7 system 
l lSouthworth et alj|201 II) . The new results are almost identical to 
the old ones (Table A125), except for the correction to pb discussed 
in Sect. 13.31 The full set of homogeneous properties for transiting 
extrasolar planetary systems can now be accessed by considering 
only the current work and Paper III. 



PERFORMANCE OF THE DEB RELATION VERSUS 
CONSTRAINTS FROM THEORETICAL MODELS 



2.5 



2.0 



One of the new procedures introduced in the current work is an 
optional constraint using a R\ = /(MA,T c ft, [tt]) relation 
obtained from we ll-studied dEBs and following the method of 
lEnoch et al.l a201Cf> . This replaces the approach used in Papers II 
and III, which utilised a mass-radius relation from dEBs, which 
was simpler but did not work very well. The cost is a reliance on 
T e ff and [|^1 measurements, which incurs a dependence on theo- 
retical model atmospheres. The new approach gives results in much 
better agreement with those found via theoretical models. 

The dEB constraint has been used to calculate physical prop- 
erties for the 30 TEPs studied in Paper III, giving a sample of 58 
TEPs with physical properties calculated in several ways: using the 
dEB constraint, using each of five different theoretical model tab- 
ulations, and the nominal results which are an unweighted mean 
of the ones from the five model sets. The parameter Kb is well- 
suited for comparing the different options, as it is the solution con- 
trol parameter in JKTABSDIM and wholly encompasses the outside 
constraints used in calculating the physical properties. A larger Kb 
results in larger numbers for all of the physical properties (see eqs. 
4 to 12 in Paper II) with the exceptions of O (which gets smaller) 
and the three quantities which have no model dependence (</b, pA 
and T c ' q ). 




0.5 



0.8 1.0 1.2 

Stellar mass (M Q ) 

Figure 38. Plot of the masses versus the radii of the stars in the 58 TEPs 
with homogeneous properties. The statistical uncertainties are shown by 
black open diamonds and the systematic uncertainties by red filled dia- 
monds. The empirical mass-radius relation from Paper II is shown with a 
blue line. 



Fig.[37] presents a detailed visualisation of the Kb values ob- 
tained from the various solutions for each TEP system. This Figure 
is a modified version of fig. 20 in Paper III with double the number 
of systems and with the dEB constraint instead of the mass-radius 
constraint. Previous assertions can be confirmed: the Claret models 
yield a larger Kb on average, the VRSS models give a lower Kb 
on average, and the Y 2 models show no trends with T c s , [^1 and 
Ma compared to the mean model solution. The dEB constraint is 
clearly hugely more successful than the mass-radius relation in re- 
producing the theoretical model results, but tends to return a larger 
Kb particularly at high metallicity. This confirms that it is a useful 
tool in quick calculation of the properties of transiting planets. 



8 PHYSICAL PROPERTIES OF THE TRANSITING 
EXTRASOLAR PLANETARY SYSTEMS 

The major results of this work are the physical properties of 32 tran- 
siting extrasolar planetary systems obtained using homogeneous 
methods and by combining all available photometric data (Tabled 
with measured spectroscopic parameters of the host stars (Table[5]l. 
The stellar properties are given in Table|6] and the planetary ones 
in Table|7] These quantities supplement (and in few cases super- 
sede) the properties for 30 objects found in Paper III, giving a total 
sample of 58 systems. The homogeneous nature of these results 
means they are well suited for comparing different TEPs, for plan- 
ning follow-up observations, and for performing detailed statistical 
studies. 

Figs. [38] and [39] show the masses and radii of the stars and 
planets with their random (black open diamonds) and systematic 
(red filled diamonds) errorbars. It is clear that the property of these 
four which is most affected by systematic error is Ma, whereas 
the masses and radii of the planets are not strongly affected by this 
model dependence, as previously found in Paper III. 

Fig.[40] shows the masses and radii of the TEPs and their host 
stars on the same plot. The plot includes all 118 TEP systems 
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Figure 37. Comparisons between the K\, values obtained using individual sets of stellar evolutionary models and the unweighted mean value, (iv"b), for each 
TEP. From left to right the panels show results for the dEB constraint and then the five stellar model sets. The top panels compare K\, to (K\,) for each model 
set, with parity indicated by a dotted line. Lower panels show the difference (K^, — (K^)) as a function of T c g , -rr and Mj±. 



known as of 201 1/04/19 and includes the 58 systems studied in this 
series of papers plus results taken from the literature for the other 
60 systems. The Sun, Jupiter, Saturn, Uranus, Neptune, Earth and 
Venus are also plotted for context, as are lines denoting the points 
where density is equal to pj up and pq. Fig.|40] clearly highlights 
the wide range of parameter space covered by these systems, as 
well as the fact that the properties of the planets are much more 
scattered than those of the parent stars. 



9 FOLLOW-UP OBSERVATIONS 

Most of the TEPs in the current work would benefit from further 
observations of some sort. In many cases the dominant uncertainty 
stems from the quality of the light curve. This remains true for 
many of the CoRoT systems, despite their space-based light curves. 
It must be remembered that CoRoT has only a 27 cm diameter 
telescope and studies relatively faint stars, so is subject to signif- 
icant photon noise. Also, several of the CoRoT TEPs have few ob- 
served transits because they were studied in short runs (CoRoT-4 
and CoRoT-15) or because they have long orbital periods (CoRoT- 



9 and CoRoT-10). Almost all of the 58 TEP systems in this series 
of papers which do not need better light curves have been observed 
from space (as well as from the ground in many cases). 

Additional RV measurements are useful too. In many cir- 
cumstances, particularly for the fainter objects, the RVs are good 
enough to unambiguously confirm the planetary nature of a system 
but are the dominant source of uncertainty in the planetary masses. 
Now over 100 TEPs are known it seems appropriate to concen- 
trate follow-up resources on measuring the physical properties of a 
golden subset of these to high precision. An additional requirement 
of RVs is definition of the orbital shape (e and w), and imprecise 
measurements of these quantities compromise measurements of the 
photometric parameters, in particular rx- 

The physical properties of quite a few of the TEPs are also 
limited by the precision of the T c g and [5^1 measurements avail- 
able. In many cases this can be improved, but in some cases this is 
not an option because the errorbars are already close to the limit set 
by our understanding of low-mass stars (taken to be 50 K in T c g 
and 0.05 dex in [2Sj ). There is no immediate prospect of lower- 
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Figure 40. Plot of the masses versus the radii of all published TEPs and their host stars. Blue crosses show the objects studied in the current series of papers 
and red crosses show results taken from the literature. The Solar system bodies are shown by green filled circles. The grey dotted lines show the loci where 
density is equal to pj up (upper line) and pQ (lower line). 
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Figure 39. Plot of the masses versus the radii of the planets in the 58 TEPs 
with homogeneous properties. The statistical uncertainties are shown by 
black open diamonds and the systematic uncertainties by red filled dia- 
monds. Blue dotted lines show where density is 1.0, 0.5 and 0.25 pj up . 



ing these threshold s; in fact there is e vidence that they are already 
slightly optimistic JBruntt e t al. 2010). 

Quite a few of the CoRoT TEPs have ephemerides which will 



become uncertain over the timescale of a few years, so need to 
be followed up soon to quash possible ephemeris drift. To inves- 
tigate this I compiled a catalogue of ephemerides of all known 
TEPs and identified the first predicted times of transit which were 
uncertain by one hour, and by half of one transit duration. A list 
of objects for which one of these dates is earlier than the year 
2020 is given in Table[9] The Kepler planets are separated because 
they continue to be observed and will have significantly improved 
ephemerides even with existing (unreleased) data, and in some 
cases have strong TTVs. The list of the other planets is clearly dom- 
inated by CoRoT objects, and it is notable that the ephemerides for 
CoRoT-4 and CoRoT-14 are already uncertain by more than one 
hour. Further photometric observations of these are advocated be- 
fore the ephemerides deteriorate much further. 



10 SUMMARY 

The physical properties of 32 transiting extrasolar planetary sys- 
tems have been derived from public light curves and published 
spectroscopic parameters of the host stars. These include 15 sys- 
tems observed by the CoRoT satellite, ten by Kepler, and five by 
the EPOCH project on the Deep Impact spacecraft. Combined with 
the 30 objects examined in Paper III, a sample of 58 TEPs with 
homogeneously measured properties is obtained. 

All available transit light curves of each TEP were obtained 
and modelled using the JKTEBOP code, with careful attention paid 
to the treatment of limb darkening, contaminating light, orbital ec- 
centricity, Poisson and correlated noise, and long effective expo- 
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Table 6. Physical properties of the stellar components of the TEPs studied in this work. For each quantity the first uncertainty is derived from a propagation of 
all observational errors and the second uncertainty is an estimate of the systematic errors arising from the dependence on stellar theory. 



System 



Semimajor axis (AU) 



Mass(M Q ) 



Radius (R Q ) 



log g A [cm/s] 



Density (p Q ) Age (Gyr) 



CoRoT-1 

CoRoT-2 

CoRoT-3 

CoRoT-4 

CoRoT-5 

CoRoT-6 

CoRoT-7 

CoRoT-8 

CoRoT-9 

CoRoT-10 

CoRoT-11 

CoRoT-12 

CoRoT-13 

CoRoT-14 

CoRoT-15 

HAT-P-4 

HAT-P-7 

HAT-P-11 

HD 17156 

HD 80606 

Kepler-4 

Kepler-5 

Kepler-6 

Kepler-7 

Kepler-8 

KOI-428 

LHS6343 

TrES-2 

TrES-3 

WASP-3 

WASP-7 

XO-4 



0.02536±0.00098±0.00016 

0.02854±0.00036±0.00032 

0.05783±0.00078±0.00035 
n nai on +0.00H0 +0.00061 

U.uyizu _ . 00112 -0.00067 

n n^nru +0-00161 +0.00022 

u.uouut _ 00092 -0.00033 
0.0855 ±0.0016 ±0.0007 
0.01690±0.00036±0.00025 
0.0633 ±0.0019 ±0.0008 
0.4027 ±0.0095 ±0.0056 
0.1060 ±0.0011 ±0.0009 
0.0440 ±0.0016 ±0.0003 
0.0394 ±0.0011 ±0.0004 
0.0510 ±0.0012 ±0.0005 

0.02687±0.00077±0.00015 

n ru^s +0.0018 +0.0005 

u - u ^ do -0.0022 -0.0003 

n nAARz. +000113 +0.00084 

u.iwtoo _ . 00062 -0.00054 

0.03805±0.00033±0.00015 
0.05259±0.00056±0.00027 
0.1637 ±0.0019 ±0.0022 
0.4564 ±0.0054 ±0.0068 
0.0449 



,+0.0024 +0.0005 
-0.0012 -0.0004 



0.04967 



0.95 ±0.11 ±0.02 

1.018±0.038±0.034 

1.403±0.056±0.026 

1 1Qzl +0.043 +0.024 
1 -0.043 -0.026 



1.025 



0.056 -0.020 
1.054±0.059±0.024 

0.884±0.056±0.039 

0.878±0.078±0.035 

0.960±0.068±0.040 

0.904±0.027±0.022 

1.26 ±0.14 ±0.02 

1.018±0.088±0.029 

1.086±0.077±0.035 

1.125±0.098±0.018 



1.31 



+0.16 +0.04 



+0.00051 +0.00014 



-0.000:5* -0.001)28 
n rizlzL^S +0.00181 +0.00080 
U.Lra^OO _ . 00081 -0.00053 

0.0613 ±0.0017 ±0.0006 
0.0485 ±0.0012 ±0.0002 
0.0795 ±0.0022 ±0.0015 
0.0850 ±0.0031 ±0.0007 
0.03567±0.00061±0.00029 
0.02276±0.00012±0.0001 1 
0.03 1 85±0.00086±0.00020 
0.0619 ±0.0010 ±0.0003 

n n^A7A +0-00162 +0.00020 
u.votlt_ 00056 -0.00031 



1.173 
1.296 



-0.19 -0.03 
1 97 -i +0.096 +0.072 
■"•••" x -0.053 -0.046 

1.511±0.039±0.017 
0.812±0.026±0.012 
1.297±0.046±0.053 
1.018±0.035±0.045 

+0.193 +0.039 
-0.095 -0.033 
+0.040 +0.011 
-0.030 -0.022 
1 114+0- 146 +0.061 
J-- 11 * -0.062 -0.040 

1.28 ±0.11 ±0.03 

1.230±0.072±0.01 

1.42 ±0.12 ±0.08 

0.440±0.049±0.012 

0.991±0.052±0.024 

0.921±0.014±0.014 

1.26 ±0.10 ±0.02 

1.285±0.063±0.019 

1 9R^+ 0117 +0-014 
1,ioJ -0.039 -0.022 



1.131 ±0.045 ±0.007 
0.907 ±0.020 ±0.010 
1.575 ±0.094 ±0.010 



4.311 ±0.019 ±0.003 
4.530 ±0.015 ±0.005 
4.191 ±0.046 ±0.003 



+0.100 +0.013 



+0.0923 +0.0077 
-0.0322 -0.0084 



1.1475 

1 nMfi+O' 0810 +0-0045 
l.uoiu _ .0666 -0.0069 

1.043 ±0.029 ±0.008 
0.96 ±0.15 ±0.01 
0.898 ±0.090 ±0.012 
0.938 ±0.059 ±0.013 
0.743 ±0.055 ±0.006 
1.374 ±0.061 ±0.009 
1.046 ±0.042 ±0.010 
1.274 ±0.077 ±0.014 
1.19 ±0.14 ±0.01 



1.36 



+0.39 +0.01 
-0.12 -0.01 



1.600 



+0.113 +0.030 
0.037 -0.019 

1.955 ±0.019 ±0.007 

0.695 ±0.014 ±0.004 

1.487 ±0.037 ±0.020 

1.037 ±0.032 ±0.015 

1 ,0+0.33 +0.02 
± - f *° -0.13 -0.01 
1 X.AA+ - 055 +0-004 
l.O'it _ .042 -0.009 
1 9R1 +0-057 +0.023 
J..ZOJ. _ .024 -0.015 

1.969 ±0.081 ±0.018 

1.495 ±0.037 ±0.005 

2.24 ±0.25 ±0.04 

0.418 ±0.030 ±0.004 

0.964 ±0.017 ±0.008 
0.8235±0.0098±0.0040 

1.366 ±0.044 ±0.008 

1.466 ±0.094 ±0.007 

+0.386 +0.006 



4.3959 



+0.0237 +0.0029 



-0.0678 -0.0032 
A An*.? +0-0683 +0.0019 
t.tvoo _ .0594 -0.0028 

4.425 ±0.022 ±0.003 
4.42 ±0.14 ±0.01 

4.475 ±0.077 ±0.006 

4.476 ±0.063 ±0.006 
4.652 ±0.062 ±0.004 
4.264 ±0.019 ±0.003 
4.407 ±0.029 ±0.004 
4.264 ±0.040 ±0.005 
4.338 ±0.082 ±0.002 



4.288 



+0.059 +0.005 
-0.191 -0.003 



+0.015 +0.008 
0.038 -0.005 

4.0354±0.0049±0.0017 

4.663 ±0.012 ±0.002 

4.207 ±0.018 ±0.006 

4.415 ±0.021 ±0.007 
+0.063 +0.005 
-0.133 -0.004 
+0.023 +0.001 
-0.024 -0.002 



4.134" 



4.168 
4.174 



4.284" 



1.531 



+0.017 +0.008 
-0.011 -0.005 

3.959 ±0.024 ±0.004 
4.178 ±0.022 ±0.002 
3.889 ±0.099 ±0.008 
4.839 ±0.035 ±0.004 
4.4660±0.0081±0.0035 
4.5710±0.0064±0.0021 
4.268 ±0.018 ±0.003 
4.215 ±0.046 ±0.002 

+0.034 +0.002 



■0.068 -0.009 



4.177" 



0.172 -0.002 



0.660 ±0.019 

1.362 ±0.064 

0.359 ± 0.058 

0.790+.°;™ 

0.88+°^ 
0.929 ±0.064 

1.00 ±0.48 
1.21 ±0.32 
1.16 ±0.24 
2.20 ±0.47 

0.488 ± 0.022 

0.889 ± 0.076 

0.526 ±0.072 

0.67 ±0.19 

52+ 012 
U - OZ -0.25 

S10+ 0016 
u.o±u_ 041 

0.2023 ±0.0024 

2.415 ±0.097 

0.395 ±0.022 

0.913 ±0.062 

362+ 0081 
u.cso^_ 136 

n ocrn + 0.025 

u,OJi -0.029 

n c w +0.0076 

u.oooo_ 0209 

0.168 ±0.012 

0.368 ±0.014 

0.126 ±0.042 

6.01 ±0.76 
1.105 ±0.011 
1.648 ± 0.041 
0.495 ± 0.024 
0.408 ± 0.068 

fl QKS+0-046 

u.ooo_ al60 



, +4.0 +0.7 
-3.8 -0.7 

hi. 9 +1.5 
-2.1 -0.6 
- +0.5 +0.3 
' -0.5 -0.2 
,+2.6 +0.6 
' -1.0 -0.3 
,+2.6 +0.6 
-5.3 -1.0 
9 c+2.1 +0.6 
z '° -1.7 -0.7 
unconstrained 

unconstrained 

unconstrained 
+2.0 +0.3 
-0.2 -0.1 

,+0.8 +0.4 
-2.1 -0.4 
,+3.3 +1.8 
' -6.7 -1.5 
,+1.4+0.5 
'-6.2 -1.0 
o y+2.5 +0.7 



7.8; 
0.6; 
1.5; 
0.8; 

3.9" 



O.i; 
2.0; 

5.8; 

5.i 



5.0 -0.6 
1 fi+ 4 - 5 + - 9 

1 -°-5.9 -1.6 
n Q +0.6 +0.6 
0,3 -0.9 -1.1 

9 n+0-4 +0.3 
z - -0.3 -0.2 
unconstrained 



2.8 ±1:1 

5.9+" 

■1.5 

-2.5 

+0.3 

-0.5 

+0.8 

-2.3 

+0.4 

-1.6 

-1.8 

3.1 

2.1+ ' 6 



5.3; 

2.8 
5.7; 
4.4; 
3.2" 



-n.o 
1+0.7 

-0.0 
1+2.0 

-2.2 

1.0+ 05 



4.0; 

3.4" 



2.1 

2.5 
2.7 



0.0 
+1.2 
-1.2 
+0.8 
-0.9 
+ 1.1 
-0.5 



+0.4 
-0.4 
+4.1 
-2.1 
+0.4 
-0.4 
+0.3 
-0.3 
+0.3 
-0.5 
+0.7 
-0.4 
+0.6 
-1.1 
+0.1 
-0.3 
+2.0 
-3.0 
+0.5 
-0.5 
+0.0 
-0.0 
+0.4 
-0.4 
+0.2 
-0.4 
+0.2 
-0.3 



sure times. The results for each light curve were then amalgamated 
to yield combined photometric parameters for the system, which 
were compared with literature results. 

The physical properties of the TEPs were calculated from 
measured quantities by applying constraints from theoretical mod- 
els, guided by the atmospheric parameters of the host stars. Five 
different sets of theoretical model tabulations were used, and the 
final results for each TEP are the unweighted mean of the individ- 
ual results for each output parameter. Systematic errors were esti- 
mated by the interagreement between the individual model results, 
and statistical errors were propagated using a perturbation analysis. 
The constants and units needed in this process were tabulated for 
reference, and an error in the unit used for planetary density was 
fixed. 

I also calculated the physical properties of each TEP system 
usin g a constraint obta ined from eclipsing binary star systems (see 
also lEnoch et alj|2010f) . The constraint was applied in the form of 
log 10 R = /(log 10 T e g, log 10 p, T^l ), where the precise equation 
and calibration coefficients were determined using the the mea- 
sured properties of 90 well-studied detached eclipsing binaries. 
This gives results in generally good agreement with those from 
using theoretical stellar models as a constraint, although a trend 
towards poorer agreement is seen at higher metallicities. It is not 



obvious whether this trend arises from an imperfection in the cal- 
ibration equation or source data, or from the physical effects in- 
cluded in the theoretical models. 

The resulting physical properties of the 32 TEP systems are 
typically in good agreement with published results, but exceptions 
exist. My results for CoRoT-5 disagree with those of the discovery 
paper, and this may be related to the treatment of orbital eccen- 
tricity. The public light curve of C0R0T-8 does not match the pub- 
lished orbital ephemeris: I measure a revised ephemeris and some- 
what different properties compared to the discovery paper. The two 
C0R0T light curves (short and long cadence) of CoRoT-13 are dis- 
crepant. After rejection of the much less reliable 512 s cadence data 
I find physical properties of the system which are very different to 
previously thought. The resulting density of the planet is almost a 
factor of three smaller, moving it from outlier status to more repre- 
sentative of the general population of transiting Hot Jupiters. Many 
of the error estimates in the literature are far smaller than I find, and 
are not supported by the intrinsic quality of the data. 

My analysis of the TEPs observed by Kepler uses data from 
Quarters 0, 1 and 2 for most of the objects. It is therefore the first 
analysis of most of the TEPs to include short-cadence data. This 
allows me to provide updated ephemerides and more reliable phys- 
ical properties. My results for Kepler-5 are somewhat different to 
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Table 7. Physical properties of the planetary components of the TEPs studied in this work. For each quantity the first uncertainty is derived from a propagation 
of all observational errors and the second uncertainty is an estimate of the systematic errors arising from the dependence on stellar theory. 



System 



Mass(M Jup ) 



Radius (Rj up ) 



9b (ms *) 



Density (p Jup ) 



T c ' q (K) 



e 



CoRoT-l 

CoRoT-2 

CoRoT-3 

CoRoT-4 

CoRoT-5 

CoRoT-6 

CoRoT-7 

CoRoT-8 

CoRoT-9 

CoRoT-10 

CoRoT-11 

CoRoT-12 

CoRoT-13 

CoRoT-14 

CoRoT-15 

HAT-P-4 

HAT-P-7 

HAT-P-11 

HD 17156 

HD 80606 

Kepler-4 

Kepler-5 

Kepler-6 

Kepler-7 

Kepler-8 

KOI-428 

LHS 6343 

TrES-2 

TrES-3 

WASP-3 

WASP-7 

XO-4 



1.03 ±0.10 ±0.01 
3.62 ±0.14 ±0.08 
21.96 ±0.65 ±0.27 



0.731 



+0.072 +0.0010 



0.073 -0.011 

+0.058 +0.004 

0.031 -0.006 

±0.34 ±0.05 



0.470 
2.96 

0.0220 ± 0.0050 ± 0.0007 

0.216 ±0.036 ±0.006 

0.826 ±0.080 ±0.023 

2.78 ±0.14 ±0.05 

2.34 ±0.39 ±0.03 

0.887 ±0.077 ±0.017 

1.312 ±0.092 ±0.028 

7.67 ±0.49 ±0.08 



64.9 



+5.3 ■ 
6.2 



1.3 
■1.0 



0.680 



+0.038 +0.026 



,+0.048 +0.006 
-0.040 -0.007 

, +0.057 +0.023 
-0.031 -0.015 



0.025 -0.016 

1.799 ±0.038 ±0.014 
0.084 ±0.007 ±0.001 
3.262 ±0.072 ±0.088 
4.114 ±0.096 ±0.122 
n ri7c;+ - 013 +0-002 

U.UIO -0.011 -0.001 

2.040 

0.633 

0.425 ±0.046 ±0.008 

0.59 ±0.12 ±0.00 

2.12 ±0.35 ±0.08 

69.9 ±5.6 ±1.2 

1.206 ±0.045 ±0.020 

1.899 ±0.060 ±0.019 

2.03 ±0.12 ±0.03 

0.96 ±0.13 ±0.01 

1 Z.A7 +0-110 +0.011 
l.O'il _ .066 -0.017 



1.551 ±0.064 

1.470 ±0.028 

1.037 ±0.069 

1 160 + 0116 
±.iou_ 041 

-I I S9 +0.102 
l.±oz _ 098 

1.185 ±0.041 
0.166 ±0.043 
0.712 ±0.083 
1.037 ±0.081 
0.941 ±0.085 
1.426 ±0.057 
1.350 ±0.074 
1.252 ±0.075 
1.018 ±0.079 

i Q45 +0.347 
J--U40 _ 091 

, oo 7 +0.075 
J-.OOI _ 32 

1.465 ±0.015 
0.397 ±0.009 
1.065 ±0.033 
1.003 ±0.023 
0.368 ±°;°£ 

-, 91n +0.035 
l.ziu _ Q 030 

1 169 + 0052 
l.iua _ 22 

1.602 ±0.075 
1.381 ±0.037 
1.27 ±0.17 
0.864 ± 0.048 
1.193 ±0.021 
1.310±0.019 
1.416 ±0.047 
1.363 ±0.093 
1.287 +g : g« 



±0.010 
±0.016 
±0.006 

+0.008 
-0.009 
+0.005 
-0.008 

±0.009 
±0.002 
±0.010 
±0.014 
± 0.008 
±0.009 
±0.013 
±0.013 
±0.005 

+0.011 
-0.008 
+0.025 
-0.016 

±0.006 
± 0.002 
±0.014 
±0.015 

+0.004 
-0.003 
+0.002 
-0.002 
+0.021 
-0.014 

±0.014 
±0.005 
±0.02 
± 0.007 
±0.010 
±0.006 
±0.009 
± 0.007 

+0.005 
-0.007 



1+240 
-620 

,+0.44 
-0.91 



10.65 ±0.69 

41.5 ±1.7 
506 ± 67 
13 5+ 16 

s.3±i:i 

52.3 ±6.4 
19 ±12 

10.6 ±2.9 
19.1 ±3.2 

78 ±14 
28.5 ±4.2 

12.1 ±1.3 

20.7 ±2.5 
183 ± 27 
1470; 
9.42; 

20.77 ±0.33 

13.2 ±1.1 
71.2 ±3.7 
101.4 ±3.9 

13.8l4g 

34.5+J> 

11.48±° ; jl 

4.10 ±0.43 

7.7 ± 1.4 

32 ± 10 

2320 ± 210 

21.02 ±0.31 

27.4 ±1.1 
25.1 ±1.2 
12.9 ±2.4 

23.ll^ 



0.259 ±0.021 ±0.002 
1.066 ±0.057 ±0.012 
18.4 ±3.7 ±0.1 

n aqq+O- 063 +0.003 
U.«fcOO _o.n7 -0.003 
fl Ofifi +0-082 +0.002 
U.ZOD _ 058 -0.001 
±0.23 



1.66 

4.5 

0.56 

0.69 

3.13 

0.76 



±0.01 

±0.1 

±0.01 

±0.01 

±0.03 

±0.00 



±4.5 

±0.21 

±0.17 

±0.8 

±0.12 
0.337 ±0.052 ±0.003 
0.62 ±0.11 ±0.01 
6.8 ±1.5 ±0.0 

c-o +13 +0 

00 -29 -0 
fl 9fifi+°- 018 +0-003 
U.ZOO_ Q3g _ .005 

0.535 ±0.011 ±0.002 
1.26 ±0.12 ±0.01 
2.52 ±0.20 ±0.03 
3.82 ±0.23 ±0.06 

-I A , +0.48 +0.01 
1-ft± -0.62 -0.02 

i n7fi+°- 080 +°- 002 

l.U/D _ 086 -0.001 

n Q7n+°- 015 +0.004 

U.J/U _ 026 -0.007 

0.097 ±0.013 ±0.001 

0.210 ±0.040 ±0.001 

0.98 ±0.44 ±0.02 

101 ±13 ±0 

0.665 ±0.015 ±0.005 

0.790 ± 0.040 ± 0.004 

0.669 ± 0.047 ± 0.004 

0.356 ±0.087 ±0.002 

n KR+ - 11 + - 00 
U - UO -0.37 -0.00 



1915 ±49 
1548 ± 22 
1695 ± 57 



1058 
1348' 



+42 
-17 
h50 

-51 

1025 ±16 
1910 ± 140 
922 ± 41 
413 ± 14 
647 ± 24 
1735 ± 34 
1410 ± 28 
1432 ± 39 
1936 ± 95 
1670" 



1+200 
-80 



imi±H 

2194 ±27 
838 ± 10 
883 ± 11 
405.0 ±7.0 
1620^° 
1692+^ 

— 25 

1451 H 



-15 

-13 

1621 ± 23 
1662 ± 41 
1666 ± 92 
352 ± 22 
1466 ± 12 
1638 ± 22 
2020 ± 35 
1502 ± 47 



0.0354 ±0.0025 ±0.0002 
0.1381 ±0.0049 ±0.0016 
1.74 ±0.12 ±0.01 
n no«9 +00099 +0.0007 

U.UJOZ _ .oi27 -0.0006 

n mss+O- 0054 +0.0003 

U.UOOO_ 0003g _ .0002 

0.405 ±0.046 ±0.003 
0.0051 ±0.0018 ±0.0001 
0.0437 ± 0.0084 ± 0.0006 
0.668 ±0.076 ±0.009 
0.693 ±0.070 ±0.006 
0.114 ±0.017 ±0.001 
0.0508 ± 0.0042 ± 0.0005 
0.0983 ±0.0080 ±0.0010 
0.360 ±0.030 ±0.002 

, n A +0Al +0.03 

*-°*_1.07 -0.05 

n n^c;7+ '00 12 +0-0004 

U.UOO( _ .0021 -0.0007 

0.0618 ±0.0010 ±0.0002 

0.0274 ±0.0022 ±0.0001 

0.772 ±0.026 ±0.010 

3.677 ±0.093 ±0.055 

+0.0025 +0.0001 
-0.0034 -0.0002 
+0.0036 +0.0002 
-0.0040 -0.0002 



0.0156 
0.1286 



0.0431 



+0.0016 +0.0005 
-0.0022 -0.0008 

0.0253 ± 0.0024 ± 0.0002 
0.0337 ±0.0063 ±0.0001 
0.188 ±0.038 ±0.003 
31.2 ±2.2 ±0.3 
0.0727 ±0.0017 ±0.0006 
0.0716 ±0.0024 ±0.0004 
0.0724 ±0.0031 ±0.0004 
0.068 ±0.010 ±0.000 

n inQq+O' 0085 +0.0006 
u.iuzo_ 00240 _ 0.0004 



those previously published, due primarily to the inclusion of the 
Quarter 2 short-cadence data. 

Asteroseismic studies are available for the three previously- 
known TEPs in the Kepler field, based on the Kepler short-cadence 
data, and for HD 17156 based on HST data. These studies use the- 
oretical stellar models to interpret the oscillation spectrum of the 
star, and measure the stellar density to very high precision. The 
corresponding values I find from the light curve analysis are in 
good agreement for HD 17156 (0.8cr) and TrES-2 (l.lcr) but not for 
HAT-P-7 (2.9a) or HAT-P-1 1 (6.5cr). This indicates a problem with 
at least one of the approaches, which might be related to underes- 
timation of the true uncertainties, starspot activity or the measured 
orbital eccentricity of the HAT-P- 1 1 system. 

Finally, the complete error budgets generated for each TEP 
system allow identification of the observations which would lead 
to the greatest improvement in our measurement of their physical 
properties. Many objects would benefit from further photometric 
observations - which continue to be obtained for the the TEPs in the 
Kepler field of view - as well as from spectroscopic radial velocity 
measurements and spectral synthesis analyses. A list of nine TEP 
systems is given whose orbital ephemerides will become uncertain 
by more than one hour within this decade; the transits of CoRoT-4 
and CoRoT-14 are already not predictable to within one hour. 

The homogeneous physical properties obtained in this work 
will be useful for detailed statistical studies of the extrasolar 



planet population as well as for planning many types of follow- 
up observations of these objects. The primary results from the 
current work and from previous papers in the series, along 
with a range of other useful information, have been concate- 
nated and placed in an online catalogue. TEPCat is available at 
|http : //www . astro . keele . ac . uk/~ jk.t/tepcat/| in a 
range of convenient formats for readers to download for reference 
and further study. 
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Table 8. Summary of which types of additional observations would be use- 
ful for the thirty TEPs studied in this work. * denotes where additional data 
would be useful, and ** indicates where it would be useful but difficult to 
either obtain or interpret. 



System 


Photometric 


Radial 


Spectral 




observations 


velocities 


synthesis 


CoRoT-1 




* 


* 


CoRoT-2 




* 


• 


CoRoT-3 


** 




* 


CoRoT-4 


** 


* 




CoRoT-5 




* 


* 


CoRoT-6 




* 


* 


CoRoT-7 


** 


** 




CoRoT-8 


* 


** 




CoRoT-9 








CoRoT-10 


** 






CoRoT-11 




** 


* 


CoRoT-12 


** 




• 


CoRoT-13 




* 




CoRoT-14 


*r* 




** 


CoRoT-15 




** 


** 


HAT-P-4 


* 






HAT-P-7 








HAT-P-11 








HD 17156 








HD 80606 








Kepler-4 


• 


* 


* 


Kepler-5 








Kepler-6 




* 




Kepler-7 


* 






Kepler-8 




* 


* 


KOI-428 


* 


* 


* 


LHS 6343 


* 


* 


• 


TrES-2 






* 


TrES-3 


* 






WASP-3 




* 


* 


WASP-7 


* 


* 


* 


XO-4 


* 
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